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Preface 



Today there exist many excellent textbooks on quantum field theory. The most popular 
ones are listed in the bibliography to the present lectures. Nevertheless, everyone who gives 
lectures on quantum field theory faces the problem of selection of material and writing the 
lecture notes for students. The present text is just the lecture notes devoted to the radiative 
corrections in QFT. On this way, one encounters two problems, namely, the ultraviolet and 
the infrared divergences. Our task is to demonstrate how one can get rid of these divergences 
and obtain finite corrections to the cross-sections of elementary processes. During the course 
we describe the methods of Feynman diagram evaluation and regularization of divergences. 
In more detail, we consider the renormalization theory and elimination of ultraviolet diver- 
gencies in the Green functions off mass shell, as exemplified by scalar and gauge theories. In 
connection with the renormalization procedure we describe also the renormalization group 
formalism in QFT. As for the infrared divergences, in the literature one can find mainly the 
discussion of the IR divergencies in quantum electrodynamics. In non-Abelian theories as 
well as in QED with massless particles the situation is much more involved as there arise 
coUinear divergences as well. In the last lecture, we show how one can get rid of these di- 
vergences using the methods developed in quantum chromodynamics. One more topic also 
related to divergences is the so-called anomalies. They also lead to unwanted ultraviolet di- 
vergent contributions. Therefore, a separate lecture is dedicated to the axial and conformal 
anomalies. 

The presented text overlaps with many textbooks and is partly borrowed from there. 
However, the composition of the material and most of the calculations belong to the author, 
so wc omit the direct references to any textbooks. It should be admitted that the style of 
presentation in different textbooks varies very much and the reader can choose the book 
according to his preferences. We mostly used the classical monograph by N.Bogoliubov and 
D.Shirkov when describing the renormalization theory and more modern book by M.Peskin 
and D.Schreder which we followed when discussing the infrared divergences. 

Our experience in giving lectures on quantum field theory, the renormalization theory and 
the renormalization group tells us that this material is still complicated for perception and is 
not always presented clearly enough. One often meets with the lack of understanding of the 
comphcated structure of the field theory which manifests itself in renormahzation theory. 
Sometimes the nonrenormalizable theories are simplistically treated as the field theories with 
a dimensional coupling constant which otherwise have no difference from the renormalizable 
ones. The collinear divergences arising in theories with massless particles, despite a long 
history, have not also become the well-known part of the QFT course. Here we make an 
attempt of a simplified presentation of this complicated material. Of course, this means that 
one has to sacrifice some rigorousness and completeness. We hope that together with the 
existing literature the present lectures will serve the goal of clarification and mastering of 
quantum field theory and its applications to particle physics. 



1 Lecture I: Radiative corrections. General analysis of 
divergent integrals 



1.1 Radiative corrections 

The formalism of quantum field theory, being the generalization of quantum mechanics to 
the case of an infinite number of degrees of freedom with nonconservation of the number of 
particles, allows one to describe the processes of scattering, annihilation, creation and decay 
of particles with the help of the set of well-defined rules. As in quantum mechanics the cross- 
section of any process is given by the square of the modulus of the probability amplitude 
calculated according to the Feynman rules for the corresponding Lagrangian integrated over 
the phase space. Since the exact calculations of the probability amplitudes seem to be 
impossible, one is bound to use the perturbation theory with a small parameter - the coupling 
constant - and get the result in the form of a power series. The leading terms of this series 
can be presented by Feynman diagrams without loops, the so-called tree diagrams. The 
examples of such diagrams for some typical processes in QED are shown in FiglH 




Figure 1: The examples of tree diagrams of different processes in QED: ) the Compton 
scattering, b) the Mueller scattering, c) the annihilation of the particle-antiparticle pair. 
Shown are the momenta of external (real) and internal (virtual) particles 

All the diagrams shown in Fig{T] are proportional to the square of the coupling constant 
e^. They are constructed according to the well-known Feynman rules and do not contain 
any integration over momenta (when working in momentum representation) since due to the 
conservation of four-momentum all momenta are defined uniquely. 

The situation changes when considering the next order of perturbation theory. As an 
example, in Figl2]we show the corresponding diagrams for the Compton scattering. 

They got the name of radiative corrections since in electrodynamics they correspond to 
the emission and absorption of photons. This name is also accepted in other theories for 
perturbative corrections. All these diagrams are proportional to the fourth power of the 
coupling constant and, hence, are the next order perturbations with respect to the tree 
diagrams. However, contrary to the tree diagrams, they contain a closed loop which requires 
the integration over the four-momenta running through the loop. Any loop corresponds to 
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Figure 2: The one- loop diagrams for the process of the Compton scattering 

the bifurcation of momenta similarly to the bifurcation of the electric current, according to 
the Kirchhoff rules, so that the total momentum is conserved but the momentum running 
along each line is arbitrary. Therefore, one has to integrate over it. 



1.2 Divergence of integrals 

Prior to calculating the radiative corrections let us consider the behaviour of the integrand 
and the integral as a whole. As an example we take the diagrams of the Compton scattering 
shown in Figl2l The integral corresponding to the diagram shown in FiglJla) has the form 



/ 



'j^{p — k + m)'-^^ 



[k'^ + ie] [{p — kY — im? + ie] ' 



where the photon propagator is written in Feynman gauge and the integration takes place in 
Minkowskian space. We shall not calculate explicitly this integral now (we shall do it later) 
but consider the integrand from the point of view of the presence of singularities as well as 
the behaviour at small and large momenta. 

The presence of poles in the propagators for momentum equal to the mass squared 
does not create any problem for the integration since according to the Feynman rules the 
denominator contains the infinitesimal imaginary term ~ e — 0, which defines the way to 
bypass the pole. The choice accepted in (11.11) corresponds to the causal Green function. 

Consider now the case of /c^ — >■ 0, the so-called infrared behaviour. Despite the presence 
of k'^ in the denominator, the singularity is absent due to the measure of the 4-dimensional 
integration which is also proportional to k'^. This is true for all such integrals. The singular- 
ities appear only for certain external momenta which are on mass shell and have a physical 
reason. Off shell the singularities are absent. For this reason we shall not discuss the infrared 
behaviour of the integrals so far. 

Consider at last the case of fc^ — * oo, the so-called ultraviolet behaviour. Notice that 
in the denominator one has 4 powers of momenta, while in the numerator one has 1 plus 
4 powers in the measure of integration. Hence one has 5-4=1, i.e. the integral is linearly 



divergent as fc^ — oo. Is it the property of a particular integral or is it a general situation? 
What happens with the other diagrams? 

Consider the integral corresponding to the diagram shown in Figj2l). One has, according 
to the Feynman rules 

f ^4, 7^(Pi - fc) + m)Y{P2 -k)+ m)-f^' 
J P[(pi-A;)2-m2][(p2-A;)2-m2]- 

We are again interested in the behaviour for oo. The counting of the powers of 

momenta in the numerator and the denominator gives: 6 in the denominator and 2 in the 
numerator plus 4 in the integration measure. Altogether one has 6-6=0, i.e., the integral is 
logarithmically divergent as fc^ — oo. 

Here we met the difficulty called the ultraviolet divergence of the integrals for the radiative 
corrections. The examples considered above are not exceptional but the usual ones. The 
corrections are infinite, which makes perturbation theory over a small parameter meaningless. 
The way out of this trouble was found with the help of the renormalization theory which will 
be considered later and now we try to estimate the divergence of the integrals in a theory 
with an arbitrary Lagrangian. 

1.3 General analysis of ultraviolet divergences 

Consider an arbitrary Feynman diagram G shown in FiglHl and try to find out whether it is 



L - the number of loops 

n_b - the number of external boson lines 

n f -the number of external fermion lines 




Figure 3: An arbitrary diagram containing L integrations 

ultraviolet divergent or not. For this purpose we have to calculate the number of powers of 
momenta in the integrand: each internal loop leads to integration d^p that gives 4 powers 
of momenta; each derivative in the vertex gives the momentum in p-space, i.e., 1; each 
internal line gives a propagator which behaves as p^^/p^, i.e., ri — 2 powers of momenta, 
where r; = 0, 1, 2 for various fields. Combining all these powers together we get the quantity 
called the index of divergence of the diagram (UV) 

a;(G)=4L+ E 5,+ (^^-2), (1-3) 

vertices internal lines 

where L is the number of loops and is the number of derivatives in a vertex v. 

The absence of the ultraviolet divergences means that uj{G) < 0. However, one has to be 
careful, there might be sub divergences in some subgraphs. Hence, the necessary condition 
for finiteness is 

The finiteness condition (UV): u^'ji) < 0, W'ji C G, 



where 7i are all possible subgraphs of the graph G including the graph G itself. 

There exists, however, a simpler way to answer the same question which does not need 
to analyse all the diagrams. One can see it directly from the form of the Lagrangian. To see 
this, let us introduce the quantity called the index of the vertex (UV) 

3 

u^ = 5^ + b., + -f^ - 4, (1.4) 

where S^, hy and are the number of derivatives, internal boson and fermion lines, respec- 
tively. Then the index of a diagram (11. 3p can be written as 

^{G)= E + 4-n,-^n^, (1.5) 

vertices 

where a;^"^ corresponds to the vertex where all the lines are internal, rib and nf are the num- 
ber of external boson and fermion lines, and we have used the fact that usually ri{hoson) = 
and ri{f ermion) = 1. 

Equation (11.51) tells us that the finiteness (a;(G) < 0) can take place if u^, < and 
the number of external lines is big enough. Prior to the formulation of conditions when it 
happens, let us consider some examples. 

Example 1: The scalar theory Cint = —X<f^- 

In this case (5„ = 0, = 0, 6^ = 4 and, hence, cj™"^' = 0. Thus, according to (II. 5p . 
uj{G) = A — rib — \'n-f and everything is defined by the number of external lines. The situation 
is illustrated in FigJH 




U)(G)=4 a)(G)=2 w(G)=0 oj(G)<0 

Figure 4: The indices of divergence of the diagrams in the scalar theory 

We see that there exists a limited number of divergent structures in the ip'^ theory. These 
are the vacuum graphs, the two- and four-point functions. All the other diagrams having 
more than 4 external lines are convergent (though may have divergent subgraphs). 

Example 2: Quantum Electrodynamics Cmt = ipAip. 

In this case 5^ = 0, = 2, = 1, ci;™"^ = 0. Hence, uj{G) = 4 — Ub — ^nj and the sit- 
uation is similar to the previous example, everything is defined by external lines. Divergent 
are the vacuum diagrams {uj{G) = 4), the photon propagator {uj{G) = 2), the electron prop- 
agator {uj{G) = 1) and the triple vertex {uj{G) = 0). All the other diagrams are convergent. 

Example 3: Four-fermion interaction Cint = Gipipipip . 



n 



Here = 0, = 4,6^ = 0, u;^"^ = 2 and, hence, uj{G) = 2N - |n/. Therefore, 
increasing the number of vertices we get new divergent diagrams independently of the number 
of external lines. The number of divergent structures happens to be infinite. 

Thus, the key role is played by the maximal index of the vertex. All the theories may be 
classified according to the value of u^'^'^ : 

{< Finite number of divergent diagrams, 
Finite number of divergent structures, (1-6) 
> Infinite number of divergent structures. 

Below we show that for the first two types of theories we can handle the ultraviolet diver- 
gences with the help of the renormalization procedure. The theories with a;™'^^ = are called 
renormalizable, the theories with cj™^^ > are called nonrenormalizable, and the theories 
with cu™'^^ < are called superrenormalizable. 



1.4 The analysis of dimensions 

The property of a theory with respect to ultraviolet divergences can be reformulated in terms 
of dimensions. Consider for this purpose an arbitrary term of the interaction Lagrangian 
which is the product of the field operators and their derivatives 



^li^) = 9Y[^iix)^^jix)■ 



'^■7) 



Consider the action which is the four-dimensional integral of the Lagrangian density 



A = J d^xC{x), 



and find the dimensions of parameters in eg. fl 1.71) . As a unit of measure we take the dimension 
of a mass equal to 1. Then the dimension of length [L] = —1, the dimension of time is also 
[T] = —1, the dimension of derivative [9^] = 1, the dimension of momenta [p^] = 1. Since 
the action is dimensionless (we use the natural units where h = c = 1) 

[A] = 0, 

the dimension of the Lagrangian is 

[C] =4, (D — in D dimensional space.) 

This gives us the dimensions of the fields. Indeed, from the kinetic term for the scalar field 
one finds 

D — 2 

[(90)^] = 4 — > [0] = 1, ( — - — in D dimensional space), 
for the spinor field 



3 
2'' 



( in D dimensional space). 



for the vector field 



,D 



in D dimensional space). 



1 n 



This allows one to find the dimension of the coupling constant in f ll.7p 



[g]=4- 6, 

Then the classification of interactions (11.61) can be written as 



;i.9) 



> Superrenormalizable theories, 
[g] = \ Renormalizable theories, 

< N onrenormalizahle thoeries. 



[1.10) 



Consider which category various theories belong to. For this purpose we have to calculate 
the dimensions of the couplings. 
Illustration 



^QED = e^lj^A^ip 



••gauge 
^Yukawa 



Thus, all these models are renormalizable. 
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[h] = -2, 
[G] = -2 

M = -1 
W = -1. 



Nonren. 
Nonren. 
Nonren. 
Nonren. 



All these models on the contrary are nonrenormalizable. Notice that they include the four- 
fermion or current- current interaction which was previously used in the theory of weak 
interactions. 

Hence, we come to the following conclusion: the only renormalizable interactions in four 
dimensions are: 

i) the ip'^ interaction; 

ii) the Yukawa interaction; 

iii) the gauge interaction; 

iv) the theory is superrenormalizable. It contains only two divergent diagrams shown 
in FigI3 





cj(G)=2 U)(G)=0 

Figure 5: The only divergent diagrams in the 0^ theory 

If one looks at the spins of particles involved in the interactions, one finds out that they 
are strongly restricted. The renormalizable interactions contain only the fields with spins 0, 



1/2 and 1. All the models with spins 3/2, 2, etc. are nonrenormalizable. The latter include 
also gravity. Indeed, the coupling constant in this case is the Newton constant which has 
dimension equal to [G] = —2, i.e., quantum gravity is nonrenormalizable. 

Since we do not know how to handle the nonrenormalizable interactions because the 
ultraviolet divergences are out of control, there are only three types of interactions which 
are used in the construction of the Standard Model of fundamental interactions, namely the 
ip^, the Yukawa and the gauge interactions with the scalar, spinor and vector particles. 

Here one has to make a comment concerning the vector fields with M ^ 0. Remind the 
form of the propagator of the massive vector field 



It gives ri = 2, which leads to some modification of the formulas used above and finally to 
the nonrenormalizability of the theory. The only known way to avoid this difficulty is the 
spontaneous breaking of symmetry. In this case. 



that gives r/ = and the theory happens to be renormalizable. This mechanism is used in 
the Standard Model to give masses to the intermediate weak bosons without breaking the 
renormalizability of the theory. 



^ M^-k^-ie ' 



■ dfif k^ky/k"^ 



1 o 



2 Lecture II: Regular izat ion 



The divergences which appear in radiative corrections are not yet a catastrophe for a theory 
(remind, for example, the infinite self-energy of an electric charge in its own Coulomb field) 
but require a quantitative description. To get a finite difference of the two infinite quanti- 
ties, one has to give them some meaning. This can be achieved by introducing a kind of 
regularization of divergent integrals. The most natural way of regularization is to cut off the 
integral on the upper or lower bound of integration. There are also different ways of regular- 
ization based on a modification of the integrand or of the measure of integration. Below we 
consider three most popular kinds of regularization: the ultraviolet cutoff in Euclidean space 
(A- regularization) , the Pauli-Villars regularization, and the dimensional regularization. 

2.1 Euclidean integral and the ultraviolet cutoff 

All the integrals in quantum field theory are written in Minkowski space; however, the 
ultraviolet divergence appears for large values of modulus of momentum and it is useful to 
regularize it in Euclidean space. Transition to Euclidean space can be achieved by replacing 
the zeroth component of momentum ko —>■ ik^, so that the squares of all momenta and 
the scalar products change the sign k"^ = k^ — k"^ ^ —kf — k"^ = —k"^ and the measure of 
integration becomes equal to d'^k — > id^kE-, where the integration over the fourth component 
of momenta goes along the imaginary axis. To go to the integration along the real axis, one 
has to perform the (Wick) rotation of the integration contour by 90° (see. Figl6]). This is 
possible since the integral over the big circle vanishes and during the transformation of the 
contour it does not cross the poles. 




Figure 6: The Wick rotation of the integration contour 



When transferring to Euclidean space the poles in all the propagators disappear. Now 
the integral in 4-dimensional Euclidean space can be evaluated in spherical coordinates and 
the integral over the modulus can be cut on the upper limit. Let us demonstrate how this 
method works in the case of the simplest scalar diagram shown in FigJTl The corresponding 

1 o 



p-k 



Figure 7: The simplest divergent diagram in a scalar theory 



pseudo-Euclidean integral has the form 



' ^ (2^ J [k^ - m^][{p - ky - m^y ^^'^^ 
Transforming it to Euclidean space one gets 

^^^^^ " (2^ / [fc| + m2][(p-fc)| + m2]- ^^-^^ 

(in what follows the index will be omitted.) 

For calculation of this kind of integrals we use the following approach. First, we transform 
the product of several brackets in the denominator into the single bracket with the help of 
the so-called Feynman parametrization. The following general formula is valid: 



1 T(ai + 0:2 H h an) 



/ dxidx2 ■ ■ ■ dxn 
Jo 



A'^'AT r(ai)r(a2) • ■ ■ FK) -^o 

5{l-Xi-X2 Xn)xT~^X2'''^ ■ ■ -^"-^ 

Here r(a) is the Euler F-function which has the following properties: 



:2.3) 



F(l) = 1, F(n + 1) = n\, xT{x) = T{x + 1), F(l + x) 



n=2 



n 



where 7^ is the Euler constant and ({n) is the Riemann zeta-function. The F-function is 
finite for positive values of the argument and has simple poles at negative integer values and 
at zero. 

In our case, {n = 2,ai = 02 = 1) and eq. fl2.3p has the form: 



1 F(2) /"i dxiX2S{l — Xi — X2 



[F + m2][(p-fc)2 + m2] F(l)F(l)yo [[k^ + m^]x^ + [{p - kf + m^]x2y 

"1 dx 



[k'^ — 2pkx + p^x + m 



212 ■ 



(2.4) 



Thus, integral 02.21) can be written as 



/f„2^ ^ ^ h j t } f d'k 

' {27i)^f J[k^-2kpx+p^x + m^y (27r)4y ""J [k2+p2x{l-x)+m^y ^ ' 



Now the integral depends only on the modulus of k and one can use the spherical coor- 
dinates: 



1 ^ /-A ^3^^ 



Hp') = 77^,idx i , , „ (2.6) 



(27r)4yo Jo +p2a;(l -x) + m2 



where the volume of the 4-dimensional sphere equals ^4 = 27r (in general — y(d/2)!- 
The integral over the modulus can be easily calculated 

1 /-A' k^dk^ _ 1 A2 

2 Jo [A;2+p2a;(l-x) +m2]2 ~ 2 °^V2a;(i _ 3,) + ^2 ) + ' ^ ■ > 

and, as one can see, is logarithmically divergent at the upper limit. The full answer has the 
form 



:2.8) 



The last integral over x can also be evaluated and takes the simplest form in the limiting 
cases for m = or p = 0. Now one can go back to Minkowski space Pe => —p^. 

The regularization with the ultraviolet cut-off is quite natural and relatively simple. The 
drawback of this regularization is Euclidean rather than Lorentzian invariance and also the 
absence of the gauge invariance. Therefore, it is not useful in the gauge theories. However, 
one should notice that the noninvariance of a regularization is acceptable since the invariance 
is restored when removing the regularization . Still, this aspect complicates the calculation 
as one has to take care of the validity of all the identities. 

2.2 Pauli-Villars Regularization 

Another method of regularization which is called the Pauli-Villars regularization is based on 
the introduction of a set of additional heavy fields with a wrong sign of the kinetic term. 
These fields are not physical and are introduced essentially with the purpose of regularization 
of divergent integrals. The main trick is in the replacement 

111 , , 

(2.9) 



p2 _ ^2 p2 _ rj^2 p2 _ ]\/[2 ' 

where M ^ cxd is the mass of the Pauli-Villars fields. As a result, the propagator for large 
momenta decreases faster, which ensures the convergence of the integrals. The divergences 
manifest themselves as logs and powers of M"^ instead of the cutoff parameter A2. 

One uses sometimes the modifications of the Pauli-Villars regularization when the re- 
placement (12. 9p is performed not for each propagator but for the loop as a whole. This 
method of regularization is called the regularization over circles. It is used in Abelian gauge 
theories for the loops made of the matter fields. This way one can preserve the gauge in- 
variance. However, in non-Abelian theories we face some problems related to the loops of 
the gauge fields which cannot become massive without violating the gauge invariance. This 
problem is often solved by introducing an additional regularization for the vector fields, for 
example, with the help of higher derivatives. Here we will not consider this regularization. 

The positive property of the Pauli-Villars regularization is the explicit Lorentz and gauge 
(in abelian case) invariance, but it requires complicated calculations since one has to calculate 
massive diagrams, while massless integrals are much simpler. 



2.3 Dimensional Regularization 



The most popular in gauge theories is the so-called dimensional regularization. In this case, 
one modifies the integration measure. 

The technique of dimensional regularization consists of analytical continuation from an 
integer to a noninteger number of dimensions. Basically one goes from some D to D — 2e, 
where e — 0. In particular, we will be interested in going from 4 to 4 — 2e dimensions. In 
this case, all the ultraviolet and infrared singularities manifest themselves as pole terms in e. 
To perform this continuation to non-integer number of dimensions, one has to define all the 
objects such as the metric, the measure of integration, the 7 matrices, the propagators, etc. 
Though this continuation is not unique, one can define a self-consistent set of rules, which 
allows one to perform the calculations. 

The metric: g^" —>■ g4^2e- Though it is rather tricky to define the metric in non-integer 
dimensions, one usually needs only one relation, namely g^^g^y = 6'^ = D = A — 2e. 

The measure: d^q — {iJ,'^yd^~'^'^q , where /i is a parameter of dimensional regularization 
with dimension of a mass. The integration with this measure is defined by an analytical 
continuation from the integer dimensions. 

The 7 matrices : The usual anticommutation relation holds {7'^,7'^} = '^^g'^"] however, 
some relations involving the dimension are modified: 

7/^7^ = D = A-2e, Tr-ff'Y = g^^'Trl = gf"" 

Usually Trl = 4 is taken. Then the 7- algebra is straightforward: 

TrYYYr = Trllg'^'^gP^ + g'^'^g'"" - g^^g^"], 

^i^^^^i^ = -7/^7^^^ + 2c/'^V = -(4 - 2e)7^ + 27^^ = -(2 - 2e)7^ etc. 

What is not well-defined is the 7^ since 7^ = 'i7°7^7^7^ and cannot be continued to an 
arbitrary dimension. This creates a problem in dimensional regularization since there is no 
consistent way of definition of 7^. 

The propagator : In momentum space the continuation is simple 

1 1 




p2 _ ^2 p2 _ j^2 



However, in coordinate space one has: (take m = for simplicity) 



p2 r^2 J p2 [x'^Y 



The basic integrals: The main idea is to calculate the integral in the space-time dimension 
where it is convergent and then analytically continue the answer to the needed dimension. 

Consider the earlier discussed example (12. ip and use the Euclidean representation (12. 5p . 
Let us rewrite it formally in D-dimensional space 



[A;2 + M2]2 2 Jo [A;2 + M2]2 



=/x(l -x) + (2.10) 



The integral over k is now the table one 

oo (kT^^-w ^^,^o_, roo ^ r(f )r(2 - f ) ^ ^2.11) 



[P + M2]2 ^ ^ Jo (X + 1)2 ' ' r(2) 

where we assume that the dimension D is such that the integral exists. In this case this is 
2 and 3. The main formula (12.111) allows one to perform the analytical continuation over 
D into the region D = 4 — 2e. For e = 0, i.e., in 4 dimensions, the integral does not exist 
since the F-function has a pole at zero argument. However, in the vicinity of zero we get a 
regularized expression. 

Collecting all together we get 

' {2n)D 2 Jo [p^x{l-x)+m^]^~D/2 ^ > 

Substituting now D = A — 2e and transforming back into the pseudo-Euclidean space one 
finds 



^ ~ (27r)4-2^ io [p2a;(l - x) ~ m^Y ^^'^^^ 
Expanding the denominator into the series over e, we finally arrive at 

= ^F(l + e) f 1 - rfa; log[^^^^i^^^] + log(47r)') . (2.14) 
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Comparing it with eg. (12.81) we see that the ultraviolet divergence now takes the form of 
the pole over e instead of the logarithm of the cutoff. This is less visual but much simpler 
in the calculations and also is automatically gauge invariant. 

We present below the main integrals needed for the one-loop calculations. They can be 
obtained via the analytical continuation from the integer values of D. We will write them 
down directly in the pseudo-Euclidean space. 
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[m? — A;2] 



The key formula is ( 12.15^ . All the rest can be obtained from it by the differentiation. Notice 
the singularity in the r.h.s. of (I2.15P for a = D/2 — n, n = 0,1, ... These integrals remain 
non-regularized. However, they usually do not appear in the real calculations. 

Let us mention one important rule used in dimensional regularization and related to the 
massless theories. By definition it is accepted that zero to any power is zero. Thus, for 
example, the following integral is zero 

d^h 

0, Va. (2.17) 



(A;2)" 
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In fact, here we have a cancellation of the ultraviolet and infrared divergences which both 
have the form of a pole over \je. There is no any inconsistency here and this way of doing 
is self-consistent in the calculations of dimensionally regularized integrals. 

This rule leads, in particular, to the vanishing of all the diagrams of the tad-pole type 
in the massless case. However, in the massive case they survive and are important for the 
restoration of the gauge invariance. As it will be clear later, in the Standard Model the 
tad-poles give their contribution to the renormalization of the quark masses and provide the 
transversality of the vector propagator in a theory with spontaneous symmetry breaking. 
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3 Lecture III: Examples of Calculations. One-loop In- 
tegrals 



All further calculations will be performed using dimensional regularization. Below we show 
how the rules described above can be applied to calculate in various models of quantum field 
theory. 



3.1 The scalar theory 

We start with the simplest scalar case and consider the theory described by the Lagrangian 

C = \{d,vf-'^v'-^^\ (3.1) 

The Feynman rules in this case are: 

i 



2 2 ' 

p — m 




-iX 



First, we find the one-loop divergent diagrams. As it follows from FigHJ they are the 
propagator of the scalar field and the quartic vertex. 

The propagator: In the first order there is only one diagram of the tad-pole type shown 
in FigE 




Figure 8: The one-loop propagator diagram 
The corresponding integral is 



Mr^)^^U^^. (3.2) 



(27r)4-2£ 2 J k'^-m? 

where 1/2 is the combinatoric factor. Calculating the integral (13.21) . according to (12.161) . we 
find 



iX Vi-l + e) 2/Ai^\F ^A 2 [l . 1 /. N , 



(3.3) 



The fact that the integral diverges quadratically manifests itself in the structure of the 
multiplier r(— 1 + e) which has a pole at e = as well as at e = 1. However, since we 
are interested in the limit e — 0, we expand the answer in the Loran series in e. As one 
can see, even in the case of quadratically divergent integrals the divergence takes the form 
of a simple pole over e, but the integral has the dimension equal to two. Notice, however. 
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that for m = the integral equals zero in accordance with the properties of dimensional 
regularization mentioned above. 

The vertex: Here one also has only one diagram but the external momenta can be adjusted 
in several ways (see FigjO]). As a result the total contribution to the vertex function consists 




Figure 9: The one-loop vertex diagram 



of three parts 

where we introduced the commonly accepted notation for the Mandelstam variables (we 
assume here that the momenta pi and p2 are incoming and the momenta and are 
outgoing) 

S= {Pl+ P2f = (PS + Pif, t = {pi- psf = {p2 - P4)^ U = {pi- Pif = {p2 - P3)^ 

and the integral equals 

' 48 (27r)4-2= J [k^ - m?][{p - kf - w?] ^ ' 

(1/48 is the combinatoric coefficient). We have already calculated this integral and the 
answer has the form (12.141) . Now we perform the calculation in a different and simpler way 
applicable to the massless integrals. 

Two comments are in order. The first one concerns the evaluation of the combinatoric 
coefficient. It comes from the expansion of the S-matrix within the Wick theorem. In the case 
when all the particles are different like, for example, in QED, the combinatoric coefficient 
is usually 1. For identical particles their permutations are taken into account already in 
the Lagrangian (the factors 1/2 and /4! in fl3.ip ) and lead to nontrivial coefficients. There 
exists a simple method to calculate the combinatoric coefficient in these cases. The coefficient 
equals 1/Sym, where Sym is the symmetry factor of a diagram. Consider the diagram shown 
in Figini If one does not distinguish the arrangement of momenta, then the diagram has 
the following symmetries: the permutation of external lines entering into the left vertex, the 
permutation of external lines entering into the right vertex, the permutation of the vertices, 
the permutation of internal lines. Altogether one has: 2x2x2x2 = 16. Hence, the 
combinatoric coefficient equals 1/16 but, since we distinguish three different momentum 
arrangements, one has 1/48. The same rule is valid for the multiloop diagrams and we will 
use it in the next section. 

The second comment is related to the calculation of the massless integrals which are much 
simpler, and in some cases one can get the answer without any explicit integration. The 
method, which we will describe below, is applicable to a certain type of massless integrals 
and is based on conformal properties of the massless integrals depending on one external 
argument and uses the symmetry between the coordinate and momentum representations. 
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The key formula is the Fourier-transformation of the propagator of a massless particle 



2 — = — ' (3-5) 

which can be generalized to an arbitrary dimension and any power of the propagator as 
follows: 

Obviously, this formula is also valid for the coordinate integration instead of momentum. 
This way the transition from momentum representation to the coordinate one and vice versa 
is performed with the help of (13.61) and is accompanied by the factor ^^^f^^^^- 

Let us go back to the diagram Figj9l In momentum space it corresponds to the integral 
over the momenta running along the loop. However, in coordinate space it is just the prod- 
uct of the two propagators and does not contain any integration. Therefore, the integral in 
momentum space can be replaced by the Fourier-transform of the square of the propagator. 
Since in the massless case all the propagators in both momentum and coordinate representa- 
tion are just the powers of p"^ or x"^, all of them are easily calculated with the help of relation 
M. 

In the case of the integral (13.41) for m = one first has to mentally transform both the 
propagators into coordinate space which, according to ( 13. 6p . gives the factor i^^^)'^, then 
multiply the obtained propagators (this gives 1 / (x^)^^^^)) and transform the obtained result 
back into momentum space that gives the factor y^^2-2£) power of momenta l/(p^)^ 

(the same as in the argument of the last F-function). Besides this, each loop contains the 
factor i{—ny^^. Collecting all together one gets 
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which coincides with (12.141) at m = 0. 

The described method for calculation of massless integrals is applicable to any integral 
depending on one external momentum (propagator type) and allows one to perform the 
calculations in any number of loops simply writing down the corresponding factors without 
explicit integration. In the case when the integral depends on more than one external 
momentum (like for a triangle or a box) and they cannot be put equal to zero the method 
is not directly applicable though some modifications are available. We do not consider them 
here. 

The four-point vertex in the one-loop approximation thus equals (we take the common 
factor l/4!</)^ out of the brackets): 



X 
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3j7.+5,og4. + i.„i^l + l.,>ii^ + l.,4 . (3.7) 



As one can see, the Euler constant and the logarithm of An always accompany the pole term 
1/e and can be absorbed into the redefinition of /i^. 



3.2 Quantum electrodynamics 

Consider now the calculation of the diagrams in the gauge theories. We start with quantum 
electrodynamics. The QED Lagrangian has the form 



QED = -\fI, + i){i-i^d^, - m)^ + e^^^A^^ - 



where the electromagnetic stress tensor is F^y = d^Ay — d^A^, and the last term in (13.81) 
fixes the gauge. In what follows we choose the Feynman or the diagonal gauge = 1). 
The Feynman rules corresponding to the Lagrangian (13. 8p are shown in FigJTOl 
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Figure 10: The Feynman rules for QED 

In quantum electrodynamics the divergences appear only in the photon propagator, the 
electron propagator, and the triple vertex. The one-loop divergent diagrams are shown in 
FigEl 




a; 



p-k p 
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Figure 11: The one- loop divergent diagrams in QED 

We begin with the vacuum polarization graph. It is given by the diagram shown in Fig. 
[TTh). The corresponding expression looks like: 



n^.(p) 



Tr[7'^(m + k)'^^{m + k — p)] 
[im? — k'^] [m^ — {k — p^] 



(3.9) 



where the "-" sign comes from the fermion loop and q = j^qfj,. We first go to dimension 
4 — 2e. Then the integral (13.91) becomes 



e^ifi'^y f ,4_2£, 7'r[7^(m + k)Y{m + k - p)] 



'^^ ^ ^(27r)4-2^y [m2-P][m2-(fc-p)2] 



(3.10) 



Let us put m = for simplicity. This will allow us to get a simple answer at the end. First, 
we calculate the trace of the 7-matrices: 



Tr-f^kYik -p)= Tr{-f''YYY)k''ik - p)" = AkP{k - pY[g^''g''" + g^" g^^ - g^^'g 
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So the integral now looks like 
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Using the Feynman parametrization and performing the integration according to the formu- 
lae given above one finds 
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/o [p2x(l-x)]= '^'^ ^'2 7o [p2x(l - j ■ 

To evaluate the remaining integrals, we use the standard integral for the Euler beta- 
function 

dxx--\l-xY-' = B{a,P) = ^^^^, 

T{2-e)T{2-e) 



which gives in our case 

C dxx^~Hl - xf~' 
Jo 

Thus, the integral (13. lip becomes 
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where we have used that r(— 1 +£:) = — y~7^- Multiplying eg. (13. 121 ) by the trace 

[g^^g""" + g^^g""^ - g'"" g'^'^Vp" = p'^p^ + p'^p'" - g^'^p^ = '^p'^p" - g^^'p^, 
[g^^g""" + g^'g""" - g^''g'"']g'"'p^ = g^^p^ + g^^p"^ - g^^i.^ - '2e)p^ = -(2 - 2e)p^g'"', 



we find 
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Expanding now over e with the help of 



Y{e) = -r(l+£), T{2-e) = (l-£)r(l-e), r(4-2e) = (3-2£)(2-2£)(l-2£)r(l-2£), 
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we finally get 
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where 
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(3.15) 



Given the expression for the vacuum polarization one can construct the photon propa- 
gator as shown in FigJT2l 




Figure 12: The photon propagator in QED 



One has 
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where n(p2) is given by eq. fl3.15p . Notice that the radiative corrections are always propor- 
tional to the transverse tensor P^,^ = g^^, — p^Pu/p^- This is a consequence of the gauge 
invariance and follows from the Ward identities. 

Consider now the electron self-energy graph FigJTTb). The corresponding integral is 
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Acting in a usual way we go to dimension 4 — 2£:, convert the indices of the 7-matrices and 
introduce the Feynman parametrization. The result is 
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(3.17) 



The integral over k can now be evaluated according to the standard formulas 

lovr"^ (47rj ^ jo [p"^x(l — — m^xj^ 

This expression can be expanded in series in e 
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(3.19) 



Notice that the linear divergence of the integral manifests itself as a simple pole in e, and 
the coefficient has the dimension equal to 1 and is Lorentz invariant (this is either p or m). 



At last, consider the vertex function FigJTTb). The corresponding integral is 
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Transfer to dimension 4: — 2e and introduce the Feynman parametrization. This gives 



rf"^(p,^) = ^!^m[dxi dy 
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The integral over k is straightforward and gives 
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2 [{p — qYy{l — x) + p2(l — x){x — y) + q'^y{x — y) — m'^xY 



As one can see, the first integral is finite and the second one is logarithmically divergent. 
Expanding in series in e we find 



rf ^-(p, q) = te-—j \^-2r- rilE - log(47r)) 
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3.3 Quantum chromodynamics 

Consider now the non-Abelian gauge theories and, in particular, QCD. The Lagrangian of 
QCD has the form 

+ d^c'^d^c^ + gr^'d^eAy, (3.24) 

where the stress tensor of the gauge field is now = d^A^ — dyA"^ + gf°-'^'^A^^A1 and the 
last terms represent the Faddeev-Popov ghosts. 

The Lagrangian (13.241) generates the following set of Feynman rules: 



the spinor propagator 
the vector propagator 
the ghost propagator 



the spino-gauge vertex 



the triple gauge vertex 



the ghost-gauge vertex 




Consider the one-loop divergent diagrams. We start with the gluon propagator. Besides 
the diagram shown in Fig JTTI) . one has additional contributions to the vacuum polariza- 
tion from the diagrams shown in FigJTSl The first diagram takes into account the gluon 
self-interaction and the second one the contribution of the Faddeev-Popov ghosts. (As has 
already been mentioned, the tad-pole diagrams should not be included since they are auto- 
matically zero.) These diagrams depend on the choice of the gauge, and to evaluate them 
we have to fix the gauge. In what follows we choose the Feynman gauge (or the diagonal 
gauge) for the gluon field. 




Figure 13: The vacuum polarization diagrams in the Yang-Mills theory 



Then for the first diagram we have the expression 
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(3.25) 



where 1/2 is a combinatorial factor and C2 is the quadratic Casimir operator which for the 
SU(N) group equals N . It comes from the contraction of the gauge group structure constants 
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Contracting the indices and going to 4 — 2e dimensions, one gets 
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(3.26) 



To calculate the integrals, one can use the formulas given above. The first step is the Feynman 
parametrization, eq. (12.41) . and then the momentum integration is performed according to 
eqs. (l2.16l) . Applying these rules we get for the integral (13.261) 
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The second diagram corresponds to the integral 
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here the "-" sign comes from the Fermi statistics of the ghost fields. 
Calculation is now straightforward and gives 
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Adding up the two contributions together, one finally has 
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or expanding in e 
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(3.31) 



Acting the same way as in QED one can calculate the contribution to the gluon propagator. 

Notice that the final result for the sum of the two diagrams is again proportional to the 
transverse tensor P^^, = g^^, — PfiPu/p"^- This is not true, however, for the diagram with the 
gauge fields and is valid only if one takes into account the ghost contribution. Notice also 
the opposite sign of the resulting expression compared to that of eq. (l3.14l) . This is due to a 
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Figure 14: The ghost propagator and the ghost-gluon vertex diagrams in QCD 



non-Abehan nature of the gauge fields and has very important consequences to be discussed 
later. 

Consider also the ghost propagator. Here there is only one diagram shown in FigJT4k). 
It corresponds to the integral 
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which equals 
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Analogously one can calculate the vertex diagrams. We consider in more detail the 
calculation of the ghost-gluon vertex as a simpler one. The corresponding diagrams are 
shown in FigJT4l To simplify the evaluation, we put one of the momenta equal to zero. Then 
the first diagram gives the integral 
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Using the equality kp = l/2[A;^+p^ — {k—pY] and substituting it into (13.341) we find that the 
first two terms are reduced to the standard integrals and the last one leads to the tad-pole 
structure and is equal to zero. Adding up all together we get 
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The second diagram gives 
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Contracting the indices in the numerator we have (p — kykp + p''k{p — k) — 2k''p{p — k), 
which after integration leads to 
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Adding up the two contributions together we find 
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Having in mind that at the tree level the vertex has the form V*"^^^ {«''c)(p) — —gfabcpp 
get the vertex function in the one-loop approximation as 
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(3.39) 
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4 Lecture IV: Renormalization. General Idea 



Thus, we have convinced ourselves that the integrals for the radiative corrections are indeed 
ultraviolet divergent in accordance with the naive power counting. The question then is: 
how to get a sensible result for the cross-sections of the scattering processes, decay widths, 
etc? To answer this question let us see what is the reason for divergences at large values 
of momenta. In coordinate space the large values of momenta correspond to the small 
distances. Hence, the ultraviolet divergences allow for the singularities at small distances. 
Indeed, the simplest divergent loop diagram (FigjT]) in coordinate space is the product of 
two propagators. Each propagator is uniquely defined in momentum as well as in coordinate 
space, but the square of the propagator has already an ill-defined Fourier-transform, it is 
ultraviolet divergent. The reason is that the square of the propagator is singular as ^ 
and behaves like In fact, the causal Green function (the propagator) is the so-called 

distribution which is defined on smooth functions. It has the (5-function like singularities 
and needs an additional definition for the product of several such functions at a single point. 
The discussed diagram is precisely this product. 

The general approach to the elimination of the ultraviolet divergences known as the TZ- 
operation was developed in the 1950s. It consists in the introduction to the initial Lagrangian 
of additional local (or quasi-local) terms, called the counter-terms, which serve the task of 
the definition of the product of distributions at the coinciding points. The counter-terms 
lead to additional diagrams which cancel the ultraviolet divergences. The peculiarity of this 
procedure, being the subject of the Bogoliubov-Parasiuk theorem, is in that the singularities 
are local in coordinate space, i.e., are the functions of a single point and can contain only 
a finite number of derivatives. In the theories belonging to the renormalizable class, where 
the number of divergent structures is finite, the number of types of the counter-terms is also 
finite, they repeat the terms of the original Lagrangian. This means that the introduction of 
the counter-terms in this case is equivalent to the modification of the coefficients of various 
terms., i.e. to the modification of the normalization of these terms. That is why this 
procedure was called the renormalization procedure. 

It should be stressed that the parameters of the original Lagrangian like the masses, the 
coupling constants and the fields themselves are not, strictly speaking, observable. They can 
be infinite. It is important that the renormalized parameters which enter the final answers 
are meaningful. 

Below we show by several examples of renormalizable theories how one introduces the 
counter-terms into the Lagrangian, how they lead to the renormalization of the original 
parameters and how the renormalization procedure allows one to get finite results for the 
Green functions. 

4.1 The scalar theory. The one- loop approximation 

We start with the one-loop approximation and consider for simplicity the scalar theory 
(13. ip . It belongs to the renormalizable type and has a finite number of ultraviolet divergent 
structures. The one-loop divergent diagrams in this theory were calculated in the third 
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lecture. Here we are interested in the singular parts, i.e., the poles in e. They are given by 
eqs. (EJl) and (O 

The propagator : Sing Ji(p^) = -^"^^( 

The vertex : Sing Ti{s,t,u) = -«A(y^)(-^). 

Note that the singular parts do not depend on momenta, i.e. their Fourier-transform has 
the form of the (5-function in coordinate space. 

In order to remove the obtained singularities we add to the Lagrangian fl3.ll) extra terms, 
the counter-terms equal to the singular parts with the opposite sign (the factor i belongs to 
the S-matrix and does not enter into the Lagrangian), namely, 

A ^ 1 A , ni^ ,o. A 3 , A ,A. , , . 

These counter-terms correspond to additional vertices shown in FigJT5| where the cross 

Figure 15: The one-loop counter-terms in the scalar theory 

denotes the contribution corresponding to (14. ip . With account taken of the new diagrams 
the expressions for the propagator (13. 3p and the vertex (13. 7p become 

W) = 3^"' (l - 7b + log(4»-) - log(mVr )) ■ (4.2) 

Notice that the obtained expressions have no infinities but contain the dependence on the 
regularization parameter /i^ which was absent in the initial theory. The appearance of this 
dependence on a dimensional parameter is inherent in any regularization and is called the 
dimensional transmutation, i.e., an appearance of a new scale in a theory. 

What we have done is equivalent to subtraction of divergences from the diagrams. In 
doing this we have subtracted just the singular parts. This way of subtraction is called the 
minimal subtraction scheme or the MS'-scheme. One can make the subtraction differently, 
for instance, subtract also the finite parts. It is useful to subtract the Euler constant and 
log47r which accompany the pole terms. This subtraction scheme is called the modified 
minimal subtraction scheme or the MS'-scheme. It is equivalent to the redefinition of the 
parameter fi"^. Another popular scheme of subtraction is the so-called MOM-scheme when 
the subtractions are made for fixed values of momenta. For example, in the case of the 
vertex function one can make the subtraction at the point s = t = u = t^. This subtraction 
is called the subtraction at a symmetric point. 

The difference between various subtraction schemes is in the finite parts; in the one- 
loop approximation this is just the constant independent of momentum, however, in higher 
loops one already has momentum dependent terms. Therefore, the finite parts of the Green 
functions depend on a subtraction scheme. Note that this dependence in general is not 



reduced to the redefinition of tfie parameter /x, since tfiere are usually a few divergent Green 
functions and all of them are independent. 

Thus, in the three subtraction schemes discussed above we have three different values for 
the vertex function 
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The counter-terms are also different. It is useful to write them in the following way 
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where for different subtraction schemes one has 
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The Lagrangian fl3.ip together with the counter-terms (14.41) can be written as 



(4.4) 
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C + AC = Z2-{d^(p) -Z — (p - Z^ — if = Csare, 



(4.6) 



where the renormalization constants Z and Z4 are given by (14. 5 p and the renormalization 
constant Z2 in the one- loop approximation equals 1. 

Writing the "bare" Lagrangian in the same form as the initial one but in terms of the 
"bare" fields and couplings 



^'b 2 -^B 4 



(4.7) 



and comparing it with (14. 6p . we get the connection between the "bare" and renormalized 
quantities 

^PB = yZ2^., m| = ZZg'^m^, \b = Z/iZ2'^\. (4.8) 



Equations (14.70 and (14.80 imply that the one-loop radiative corrections calculated from the 
Lagrangian (14.70 with parameters chosen according to (14.8114.50 are finite. 
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4.2 The scalar theory. The two-loop approximation 



Consider now the two-loop diagrams. For simplicity and in order to complete all the integra- 
tions we restrict ourselves to the massless case. Since we are going to calculate the diagrams 
off mass shell, no infrared divergences may appear. 

The propagator: In this order of PT there is only one diagram shown in FigJT6l 




Figure 16: The two-loop propagator type diagram 
The corresponding integral equals 
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(1/3! is a combinatorial coefficient). Let us use the method of evaluation of the massless 
diagrams described above. One has to transform each of the propagators into coordinate 
space, multiply them and transform back to momentum space. This reduces to writing down 
the corresponding transformation factors. One gets 
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where the Euler constant and log47r are omitted. 

The appeared ultraviolet divergence, the pole in e, can be removed via the introduction 
of the (quasi)local counter-term 



A£ = l(Z2-l)(a# 



(4.9) 



where the wave function renormalization constant Z2 in the MS scheme is obtained by taking 
the singular part of the integral with the opposite sign 
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After that the propagator in the massless case takes the form 



(4.10) 
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+ crossed terms 



Figure 17: The two- loop vertex diagrams 

The vertex: In the given order there are two diagrams (remind that in the massless case 
the tad-poles equal to zero) shown in FigJTTl 

The first diagram by analogy with the one- loop case equals the sum of s,t and u channels 

hi = /2l(s) +/2l(t) +/2l(n), 

where each integral is nothing else but the square of the one-loop integral 

^''^'^ - [j2^' J k^p-ky ) - "96 (Te^^i + 2 + In — ) . (4.12) 

(1/96 is the combinatorial coefficient). 

Opening the bracket we, for the first time here, come across the second order pole term 
1/e^ and the single pole \og{— fi"^ / s) / e accompanying it. This latter pole is not harmless 
since its Fourier-transform is not a local function of coordinates. This means that it can not 
be eliminated by a local counter-term. This would be an unremovable problem if it were 
not the one-loop counter-terms (14.11) which created the new vertices shown in FigJTSl In the 
same order of one gets additional diagrams presented in FigJTSl 




Figure 18: The diagrams with the counter-terms in the two-loop approximation 

These diagrams lead to the subtraction of divergences in the subgraphs (left and right) 
in the first diagram of FigJTTl The subtraction of divergent subgraphs (the 7?--operation 
without the last subtraction called the 7?.'-operation) looks like 




where the subgraph surrounded with the dashed line means its singular part, and the rest 
of the graph is obtained by shrinking down the singular subgraph to a point. The result has 
the form 

7^'/2l(s) = --7-r^ - + 2 + ln^ 2 + ln^ = 

4 (167r"^)^ ye —s e e ) 



Notice that after the subtractions of subgraphs the singular part is local, i.e. in momentum 
space does not contain \np^. The terms with the single pole 1/e are absent since the diagram 
can be factorized into two diagrams of the lower order. 

The contribution of a given diagram to the vertex function equals 

+ 1„2 i^l + 4 In + W + 4 In + In' i^l + 4 111 i^lH 

—s —s —t —t —u —u j 



The contribution to the renormalization constant of the four-point vertex in the MS scheme 
is equal to the singular part with the opposite sign 

The second diagram with the crossed terms contains 6 different cases. Consider one of 
them. Since we are interested here in the singular parts contributing to the renormalization 
constants, we perform some simplification of the original integral. We use a very important 
property of the minimal subtraction scheme that the renormalization constants depend only 
on dimensionless coupling constants and do not depend on the masses and the choice of 
external momenta. Therefore, we put all the masses equal to zero, and to avoid artificial 
infrared divergences, we also put equal to zero one of the external momenta. Then the 
diagram becomes the propagator type one: 



The corresponding integral is: 
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48 (27r)8-4^ J q'^(k-q)^k^(p-k) 

(1/48 is the combinatorial coefficient). Since putting one of the momenta equal to zero we 
reduced the diagram to the propagator type, we can again use the advocated method to 
calculate the massless integral. One has 
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As one can see, in this case we again have the second order pole in e and, accordingly, 
the single pole with the logarithm of momentum. The reason of their appearance is the 
presence of the divergent subgraph. Here we again have to look at the counter-terms of the 
previous order which eliminate the divergence from the one-loop subgraph. The subtraction 
of divergent subgraphs (the 7^-operation without the last subtraction) looks like 
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Once again, after the subtraction of the divergent subgraph the singular part is local, i.e. in 
momentum space does not depend on lnp2. 

The contribution to the vertex function from this diagram is: 



and, accordingly. 
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Thus, due to (14. 5 p and (I4.16P in the two-loop approximation the quartic vertex renor- 
malization constant in the MS scheme looks like: 



With taking account of the two-loop renormalization of the propagator (]4.10p one has: 

Z, = Z,Z^' = 1 + 1^ + ( ( A - — )• (4-18) 

2 2£l67r2 \167r2; M£2 12£^ ^ ' 

The statement is that the counter-terms introduced this way eliminate all the ultraviolet 
divergences up to two-loop order and make the Green functions and hence the radiative 
corrections finite. In the case of nonzero mass, one should also add the mass counter-term. 



4.3 The general structure of the R-operation 



We are ready to formulate now the general procedure of getting finite expressions for the 
Green functions off mass shell in an arbitrary local quantum field theory. It consists of: 

In any order of perturbation theory in the coupling constant one introduces to the La- 
grangian the (quasi) local counter-terms. They perform the subtraction of divergences in the 
diagrams of a given order. The subtraction of divergences in the subgraphs is provided by the 
counter-terms of the lower order. After the subtraction of divergences in the subgraphs the 
rest of the divergences are always local. The Green functions of the given order calculated 
on the basis of the initial Lagrangian with account of the counter-terms are ultraviolet finite. 

The structure of the counter-terms as functions of the field operators depends on the 
type of a theory. According to the classification discussed in the first lecture, the theories 
are divided into three classes: superrenormalizable (a finite number of divergent diagrams), 
renormalizable (a finite number of types of divergent diagrams) and non-renormalizable (a 
infinite number of types of divergent diagrams). Accordingly, in the first case one has a 
finite number of counter-terms; in the second case, a infinite number of counter-terms but 
they repeat the structure of the initial Lagrangian, and in the last case, one has an infinite 
number of structures with an increasing number of the fields and derivatives. 

In the case of renormalizable and superrenormalizable theories, since the counter-terms 
repeat the structure of the initial Lagrangian, the result of the introduction of counter-terms 
can be represented as 

£ + A£ = Cnare = ^^B, {qb}, W}), (4.19) 

i.e., Csare is the same Lagrangian C but with the fields, masses and coupling constants being 
the "bare" ones related to the renormalized quantities by the multiplicative equalities 

0f- = Z^/2((^},l/£)0, gf<^re^zii{g},l/e)g„ mf-^ = l/£)m„ (4.20) 

where the renormalization constants Zj depend on the renormalized parameters and the 
parameter of regularization (for definiteness we have chosen 1/e). In some cases the renor- 
malization can be nondiagonal and the renormalization constants become matrices. 

The renormalization constants are not unique and depend on the renormalization scheme. 
This arbitrariness, however, does not infiuence the observables expressed through the renor- 
malized quantities. We will come back to this problem later when discussing the group of 
renormalization. In the gauge theories Zi may depend on the choice of the gauge though in 
the minimal subtraction scheme the renormalizations of the masses and the couplings are 
gauge invariant. 

In the minimal schemes the renormalization constants do not depend on dimensional 
parameters like masses and do not depend on the arrangement of external momenta in the 
diagrams. This property allows one to simplify the calculation of the counter-terms putting 
the masses and some external momenta to zero, as it was exemplified above by calculation 
of the two-loop diagrams. In making this trick, however, one has to be careful not to create 
artificial infrared divergences. Since in dimensional regularization they also have the form 
of poles in s, this may lead to the wrong answers. 

In renormalizable theory the finite Green function is obtained from the "bare" one, i.e., 
is calculated from the "bare" Lagrangian by multiplication on the corresponding renormal- 
ization constant 

r({p'}, /^', g^.) = Zr{l/e, y^)rBare({p'}, gBare): (4.21) 



where in the n-th order of perturbation theory the "bare" parameters in the r.h.s. have to 
be expressed in terms of the renormahzed ones with the help of relations (14.201) taken in the 
(n-l)-th order. The remaining constant Zr creates the counter-term of the n-th order of the 
form A£ = {Zr — l)Or, where the operator Or reflects the corresponding Green function. If 
the Green function is finite by itself (for instance, has many legs) , then one has to remove the 
divergences only in the subgraphs and the corresponding renormalization constant = 1. 

Note that since the propagator is inverse to the operator quadratic in fields in the La- 
grangian, the renormalization of the propagator is also inverse to the renormalization of the 
1-particle irreducible two-point Green function 

The propagator renormalization constant is also the renormalization constant of the corre- 
sponding field, but the fields themselves, contrary to the masses and couplings, do not enter 
into the expressions for observables. 

We would like to stress once more that the 7^-operation works independently on the 
fact renormalizable or non-renormalizable the theory is. In local theory the counter-terms 
are local anyway. But only in renormalizable theory the counter-terms are reduced to the 
multiplicative renormalization of the finite number of fields and parameters. 

One can perform the 7^-operation for each diagram separately. For this purpose one has 
first of all to subtract the divergences in the subgraphs and then subtract the divergence 
in the diagram itself which has to be local. This serves as a good test that the divergences 
in the subgraphs are subtracted correctly. In this case the 7?--operation can be symbolically 
written in a factorized form 

nG= n (l-M^)G, (4.23) 

div. subgraphs 

where G is the initial diagram, M is the subtraction operator (for instance, subtraction of the 
singular part of the regularized diagram) and the product goes over all divergent subgraphs 
including the diagram itself. By a subgraph we mean here the 1-particle irreducible diagram 
consisting of the vertices and lines of the diagram which is UV divergent. The 1-particle 
irreducible is called the diagram which can not be made disconnected by deleting of one line. 

We have demonstrated above the application of the 7?.-operation to the two-loop diagrams 
in a scalar theory. Consider some other examples of diagrams with larger number of loops 
shown in FigJTOl They appear in the 0*^ theory in the three-loop approximation. 
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Figure 19: The multiloop diagrams in the 0^ theory 

In order to perform the 7^-operation for these diagrams one first has to find out the 
divergent subgraphs. They are shown in Figj20l 
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Figure 20: The divergent subgraphs in the diagrams of FigJTOl 



Let us use the factorized representation of the 7^-operation in the form of f l4.23p . For the 
three chosen diagrams one has, respectively, 
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where 71 and 72 are the one- and two-loop divergent subgraphs shown in Figj20l 

The result of the application of the 7^-operation without the last subtraction ( TZ'- 
operation) for the diagrams of interest graphically is as follows: 




«' xO< = xO< - {O: x:>cx -<Q:^- ::0>) )0< 

Figure 21: The 7?.'-operation for the multiloop diagrams 

Here, as before, the graph surrounded with the dashed circle means its singular part and 
the remaining graph is obtained by shrinking the singular subgraph to a point. 

Let us demonstrate how the T^'-operation works for the diagram FigJT^). Since the 
result of the T^'-operation does not depend on external momenta, we put two momenta on 
the diagonal to be equal to zero so that the integral takes the propagator form. Then we 
can use the method based on Fourier-transform, as it was explained above. One has 
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We use here the angular integration measure in the 4 — 2e dimensional space accepted 
above, which results in the multiplication of the standard expression by r(l — e) in order to 
avoid the unwanted transcendental functions. Following the scheme shown in Figj2T] we get 
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Combining all together one finds 




1 

£»(l-2e) Vp2 ; 



Note the cancellation of all nonlocal contributions. The singular part after the i?'-operation 
is always local. 

The realization of the T^'-operation for each diagram G allows one to find the contribution 
of a given diagram to the corresponding counter-term and, in the case of a renormalizable 
theory, to find the renormalization constant equal to 

Z = 1-IC n'G, (4.24) 

where /C means the extraction of the singular part. Adding the contribution of various 
diagrams we get the resulting counter-term of a given order and, accordingly, the renormal- 
ization constant. 
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5 Lecture V: Renormalization. Gauge Theories and 
the Standard Model 

Consider now the gauge theories. The difference from the scalar case is in the relations 
between various renormalization constants which follow from the gauge invariance. If the 
regularization and the renormalization scheme do not break the symmetry these relations 
hold automatically. In the opposite case, this is an additional requirement imposed on the 
counter-terms. 



5.1 Quantum electrodynamics 

Quantum electrodynamics (13. 8p is a renormalizable theory; hence, the counter-terms repeat 
the structure of the Lagrangian. They can be written as 

^C-QED = -^^i^i + (^2 - 1)^^^^ - m{Z - 1)^^ + e(Zi - (5.1) 

The term that fixes the gauge is not renormalized. In the leading order of perturbation 
theory we calculated the corresponding diagrams with the help of dimensional regulariza- 
tion (see (13. 15p . (13. 191) . (13. 231) ). Their singular parts with the opposite sign give the proper 
renormalization constants. They are, respectively, 

Z - I ^ 
^1 — i-- 77—^-, 
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Adding dSH]) with (^M> we get 



Cqed + ^Cqed = -^F^u + Z2itpdtP - mZtPi; + eZi^i^ - ^{d^A^f 

'\fIuB + ii^sdipB - mZZ^^ijBipB + eZiZ^^Z^^^^tpBABtpB 
Z^' 



id,A^B)\ (5.3) 



2e 

that gives 

ipB = Zl^'^i/j, Ab = Zlj'^A, VTLB = ZZ^^m, cb = ZiZ^^ Z^^^'^e, ^b = Z^^. (5.4) 



The gauge invariance here manifests itself in two places. First, the transversality of the 
radiative correction to the photon propagator means that the gauge fixing term is not renor- 
malized and, hence, the gauge parameter ^ is renormalized as a gauge field. Second, the 
gauge invariance connects the vertex Green function and the fermion propagator (the Ward 
identity), which leads to the identity Z\ = Z^. Since the dimensional regularization which 
we use throughout the calculations does not break the gauge invariance, this identity is sat- 
isfied automatically (see 05.21) ). This means that the renormalization of the coupling 05.41) 
is defined by the photon propagator only. Note, however, that this is not true in general in 
a non-Abelian theory. 

5.2 Quantum chromodynamics 

The complications which appear in non-Abelian theories are caused by the presence of many 
vertices with the same coupling as it follows from the gauge invariance. Hence, they have 
to renormalize the same way, i.e there appear new identities, called the Slavnov- Taylor 
identities. The full set of the counter-terms in QCD looks like 

A^Z^-\)^^r^'^r''^AlAlAlAl^{Z^-\)d^c^^^^ 

+t{Z2 - l)tpdtp - m{Z - V)^^ + g{Z^^ - V)i)k'T''^, (5.5) 
that being added to the initial Lagrangian gives 

Lqj, + ^Lqb = -^{d,A: - d^A^f - gZ,r'^AlAld,Al 
-zJ^r'^^r^^AlAlAlAl - Z.d.edd' - gz.r'^d.eAy 
+tZ2^dij - mZt/j^ + (7Ziv,^i"T> - ^{d^A''^)^ 
= -^id^At^ - d^A^^r - gZ,Zf'^r'^Alj,AUd,Al^ 

-Z,Zf^^r'^r''^A'^j,Alj,A%Alj, + d,-c%d,c% + gZ,Z^'Z,''^r'^d,-c%A'^j,c% 

+ ^{d,Alj,f + tiJsd^lJB - mZZ:,^^s^B + gZ^.^Z:,^Z^^I^BhT'''^B. (5.6) 

This results in the relations between the renormalized and the "bare" fields and couplings 

V'B = Z]!"^^, Ab = Zl^'^A, cb = Zl^'^c, 

me = ZZ^'^m, ge = ZiZ^^^^g, ^b = Z^^, (5.7) 
Z\ Z^ = Z\ Z^ , Z4. = Z-^ Z^ , Z\^ ^2 ~ Z\ z^ 

The last line of equalities follows from the requirement of identical renormalization of the 
coupling in various vertices and represents the Slavnov- Taylor identities for the singular 
parts. 
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The explicit form of the renormahzation constants in the lowest approximation follows 
from the one- loop diagrams calculated earlier (see (13.141) . (I3.19p . (I3.23p . (13.311) . (I3.33p . (I3.39p . 
Aa usual, one has to take the singular part with the opposite sign. One has in the MS scheme 
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Zi = I -4-.% (5-8) 
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where the following notation for the Casimir operators of the gauge group is used 

jabcjdbc ^ fj^^ad^ {T''T'')ij = CpSij, TriT^T^) = Tpd^K 

For the SU{N) group and the fundamental representation of the fermion fields they are 
equal to 

-1 1 
Ca = N, Cf = ^^, Tf = -. 



5.3 The Standard Model of fundamental interactions 

In the Standard Model of fundamental interactions besides the gauge interactions and the 
quartic interaction of the Higgs fields there are also Yukawa type interactions of the fermion 
fields with the Higgs field. These interactions are also renormalizable and is characterized 
by the Yukawa coupling constants, one for each fermion field. The peculiarity of the SM is 
that the masses of the fields appear as a result of spontaneous symmetry breaking when the 
Higgs field develops a vacuum expectation value. As a result the masses are not independent 
but are expressed via the coupling constant multiplied by the vacuum expectation value. 
Here there are two possibilities: to treat the Yukawa couplings as independent quantities 
and to renormalize them in a usual way and then express the renormalized masses via the 
renormalized couplings or to start with the masses of particles and to treat the Yukawa 
couplings as secondary quantities. The first approach is usually used within the minimal 
subtraction scheme where the renormalizations do not depend on masses. On the contrary, 
in the scheme when the subtraction is carried out on mass shell (the so-called "on-shell" 
scheme), one usually takes masses of particles as the basis. Under this way of subtraction 
the pole of the propagator is not shifted and the renormalized mass coincides with the mass 
of a physical particle. Below we consider the renormalizations in the SM in the MS scheme 
and concentrate on the renormahzation of the fields and the couplings. 

Another property of the Standard Model is that it has the gauge group SUc{S) x SUl(2) x 
[/y (1) which is spontaneously broken to S'[/c(3) x U em(1)- In the theories with spontaneously 
broken symmetry, according to the Goldstone theorem there are massless particles, the 
goldstone bosons. These particles indeed are present in the SM but they are not the physical 
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degrees of freedom and due to the Higgs effect are absorbed by vector bosons turning into 
longitudinal degrees of freedom of massive vector particles. 

Thus, there are two possibilities to formulate the SM as a theory with spontaneous 
symmetry breaking: the unitary formulation in which nonphysical degrees of freedom are 
absent and vector bosons have three degrees of freedom, and the so-called renormalizable 
formulation in which goldstone bosons are present in the spectrum and vector fields have 
two degrees of freedom. These two formulations correspond to two different choices of the 
gauge in spontaneously broken theory. 

In unitary gauge we have only physical degrees of freedom, i.e., the theory is automatically 
unitary, hence the name of this gauge. However, the propagator of the massive vector fields 
in this case has the form 

i.e., does nor decrease when momentum goes to infinity. This leads to the increase in the 
power of divergences and the theory happens to be formally nonrenormalizable despite the 
coupling constant being dimensionless. We have mentioned this fact in the first lecture. 

On the other hand, in renormalizable gauge, where the vector fields have two degrees of 
freedom, the propagator behaves as 

which obviously leads to a renormalizable theory which explains the name of this gauge. 
However, the presence of the goldstone bosons calls into question the unitarity of the theory 
since transitions between the physical and unphysical states become possible. 

Since all the gauges are equivalent, one can work in any of them but in the unitary 
gauge one has to prove the renormalizability while in the renormalizable gauge one has to 
prove unitarity. The gauge invariance of observables preserved in a spontaneously broken 
theory should guarantee the fulfilment of both the requirements simultaneously. Note that 
in spontaneous symmetry breaking the symmetry of the Lagrangian is preserved, it is the 
boundary condition that breaks the symmetry. 

The rigorous proof of that the theory is simultaneously renormalizable and unitary is 
not so obvious and eventually was awarded the Nobel prize, but can be seen by using some 
intermediate gauge called the i?^-gauge. The gauge fixing term in this case is chosen in the 
form 

where v is the vacuum expectation value of the Higs field, and Xi the goldstone bosons. 
In this gauge the vector propagator has the form 

and at ^ = corresponds to the renormalizable gauge while as ^ — > oo it corresponds to 
the unitary one. Since all the observables do not depend on ^, we can choose ^ = when 
investigating the renormalizability properties and choose ^ = oo in examining the unitarity. 
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Since we are interested here in the renormahzabihty of the SM, in what follows we will work 
in a renormalizable gauge. 

The Lagrangian of the Standard Model consists of the following three parts: 

^ ^ gauge ~l~ ^Yukawa ~l~ ^Higgsi (^■^) 

The gauge part is totally fixed by the requirement of the gauge invariance leaving only the 
values of the couplings as free parameters 

^gauge = -^G^^G"^^ " -^Iv^lv " -^B^^B^^ (5.10) 

where the following notation for the covariant derivatives is used 

Ga a r^a q ria , ^ rabc/~ib r~ic 

iiv — (Jfi'^v ~ '^ly'^fi + 9sJ ^f^^u^ 

= d,Wl-d,W; + ge^^''W/,Wl:, 
— d^Bi, — di^B/j,, 

D^E^ = {d^ + ig'B^)Ea, 

D,Qa = {d.-i^-T'Wl-i^-B.-i^-^X^GDQ^, 

D^Ua = {d.-i^B^-i^X'^GDU^, 

D^D^ = {d^ + i^-g'B^-i^X^Gl)D^. 

The Yukawa part of the Lagrangian which is needed for the generation of the quark 
and lepton masses is also chosen in the gauge invariant form and contains arbitrary Yukawa 
couplings (we ignore the neutrino masses, for simplicity) 

^Yukawa = VafiLaEpH + v^pQ^DpH + y'^pQJJpH + h.c, (5.11) 
where H = iT2H^ ■ 

At last the Higgs part of the Lagrangian contains the Higgs potential which is chosen 
in such a way that the Higgs field acquires the vacuum expectation value and the potential 
itself is stable 

Cmggs ^-V^ m'H^H - ^{H'^Hf. (5.12) 

Here there are two arbitrary parameters: A. The ghost fields and the gauge fixing terms 
are omitted. 

The Lagrangian of the SM contains the following set of free parameters: 

• 3 gauge couplings g^, g, g'; 

• 3 Yukawa matrices y^^, y^^, y^^; 

• Higgs coupling constant A; 



• Higgs mass parameter m^; 

• the number of the matter fields (generations). 

All particles obtain their masses due to spontaneous breaking of the SUieft{^2) symmetry 
group via a nonzero vacuum expectation value (v.e.v.) of the Higgs field 

<if>=|^Q^, v = m/VX. (5.13) 

As a result, the gauge group of the SM is spontaneously broken down to 

SUc{3) ® SUl{2) ® f/y(l) SUc{3) ® Uem{1). 
The physical weak intermediate bosons are linear combinations of the gauge ones 

W^ = ^^-j=-^, Z^ = -smewB^ + cosewW^ (5.14) 

with masses ^ 

mw = —/=gv, mz = mw/ cos 9w, ^a.n9w = g'/g, (5.15) 
V 2 

while the photon field 

7/, = cosOwB^ + sinOwW^ (5.16) 

remains massless. 

The matter fields acquire masses proportional to the corresponding Yukawa couplings: 
M^p = y^pv, = yipv, = ylpv, mj, = V2m. (5.17) 

The mass matrices have to be diagonalized to get the quark and lepton masses. 

The explicit mass terms in the Lagrangian are forbidden because they are not SUieft{'2) 
symmetric. They would destroy the gauge invariance and, hence, the renormalizability of the 
Standard Model. To preserve the gauge invariance we use the mechanism of spontaneous 
symmetry breaking which, as was explained above, allows one to get the renormalizable 
theory with massive fields. 

The Feynman rules in the SM include the ones for QED and QCD with additional new 
vertices corresponding to the SU (2) group and the Yukawa interaction, as well as the vertices 
with goldstone particles if one works in the renormalizable gauge. We will not write them 
down due to their complexity, though the general form is obvious. 

Consider the one-loop divergent diagrams in the SM. Besides the familiar diagrams in 
QED and QCD discussed above one has the diagrams presented in FigJ2^ The diagrams 
containing the goldstone bosons are omitted. The calculation of these diagrams is similar 
to what we have done above. Therefore, we show only the results for the renormalization 
constants of the fields and the coupling constants. They have the form (for the gauge fields 
we use the Feynman gauge) 
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Figure 22: Some divergent one-loop diagrams in the SM. The dotted hne denotes the Higgs 
field, the solid line - the quark and lepton fields, and the wavy line - the gauge fields 
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where, for simplicity, we ignored the mixing between the generations and assumed the 



A T 



Yukawa matrices to be diagonal. 

The difference from the expressions considered above is that the renormahzation constant 
of the scalar coupling contains the terms of the type g'^/X and This is because writing 

the counter-term for the quartic vertex we factorized A. The counter-terms themselves are 
proportional to g'^ and and are not equal to zero. Thus, the quantum corrections generate 
a new interaction even if it is absent initially. Since the gauge and Yukawa interactions 
belong to the renormalizable type, the number of types of the counter-terms is finite and the 
only new interaction which is generated this way, if it was absent, is the quartic scalar one. 
With allowance for this interaction the model is renormalizable. 

Since the masses of all the particles are equal to the product of the gauge or Yukawa 
couplings and the vacuum expectation value of the Higgs field, in the minimal subtraction 
scheme the mass ratios are renormalized the same way as the ratio of couplings. To find the 
renormahzation of the mass itself, one should know how the v.e.v. is renormalized or find 
explicitly the mass counter-term from Feynman diagrams. In this case, one has also to take 
into account the tad-pole diagrams shown in Figj22l including the diagrams with goldstone 
bosons. 

For illustration we present the renormahzation constant of the 6-quark mass in the SM 



The result for the t-quark can be obtained by replacing h by t. For the light quarks the 
Yukawa constants are very small and can be ignored in eq. fl5.18l) . 

Note that here we again have the Higgs self-interaction coupling A in the denominator. It 
appears from the tad-pole diagrams but, contrary to the previous case, the renormahzation 
constant Zm^ is not multiplied by A and the denominator is not cancelled. This does not 
lead to any problems in perturbation theory since by order of magnitude A ~ (7^ ~ j/^ and 
the loop expansion is still valid. 
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6 Lecture VI: Renormalization Group 

The procedure formulated above allows one to eliminate the ultraviolet divergences and get 
the finite expression for any Green function in any local quantum field theory. In renormal- 
izable theories this procedure is reduced to the multiplicative renormalization of parameters 
(masses and couplings) and multiplication of the Green function by its own renormalization 
constant. This is true for any regularization and subtraction scheme. Thus, for example, in 
the cutoff regularization and dimensional regularization the relation between the "bare" and 
renormalized Green functions looks like 

T{{p'],^i\{g^}) = Zr{^V^^\{9^.WBare{{p'],^A9Bare]) (6.1) 
r({p'},/i', {9,}) = Zv{l/e. {9>.WBare{{p'}. V^, {9Bare}), (6.2) 

where {p^} is the set of external momenta, {g} is the set of masses and couplings, and 

gSare = Z g{{h? / , {9tJi})9 9Bare = Zg{{l/e, {g^,})g. 

It is obvious that the operation of multiplication by the constant Z obeys the group 
property. Indeed, after the elimination of divergences one can multiply the couplings, masses 
and the Green functions by finite constants and this will be equivalent to the choice of 
another renormalization scheme. Since these finite constants can be changed continuously, 
we have a continuous Lie group which got the name of renormalization group. The group 
transformations of multiplication of the couplings and the Green functions are called the 
Dyson transformations. 

6.1 The group equations and solutions via the method of charac- 
teristics 

In what follows we stick to dimensional regularization and rewrite relation (16.21) in the form 

TBareiip'}, l/s, {gBare}) = Zr\l/e, {(?,.})!({/}, f,^ , {g^}). (6.3) 

It is obvious that the l.h.s. of this equation does not depend on the parameter of dimensional 
transmutation /i and, hence, the r.h.s. should not also depend on it. This allows us to write 
the functional equation for the renormalized Green function. Differentiating it with respect 
to the continuous parameter /i one can get the differential equation which has a practical 
value: solving this equation one can get the improved expression for the Green function 
which corresponds to summation of an infinite series of Feynman diagrams. 

Consider an arbitrary Green function F obeying equation (16. 2p with the normalization 
condition 

F({/},;,2^0) = 1. 

Differentiating (16.21) with respect to /i^ one gets: 

2d ( 2 ^ 2 ^9 ^ \ -r 2d^^Zr 

^dir^^ = [^dj^^^ W^Wg) ^ = ^ -di^^^^^-^- 
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or 

(^'a^ + m§-g + 7r) r({p2}, g,) = 0, (6.4) 

where we have introduced the so-called beta function (3{g) and the anomaly dimension of 
the Green function 7r(5') defined as 



.2 dg 
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7r = (6.6) 

Equation fl6.4p is called the renormalization group equation in partial derivatives (in Ovsyan- 
nikov form). In the western literature it is also called the Callan-Simanzik equation. 

The solution of the renormalization group equation can be written in terms of character- 
istics: 

t 

r(e'M<,).r(M,,(i,9))ei , (6.7) 

where the characteristic equation is (for definiteness we restrict ourselves to a single coupling) 

j^git,g)=m, 9(0,9) = 9- (6-8) 

The quantity g(t,g) is called the effective charge or effective coupling. 

We will consider the useful properties of this solution (16. 7p later and we first derive several 
other similar equations. Since the vertex function usually comes with the coupling, one can 
consider the product 

9^ [^^9). (6.9) 
If r is the n-point function, then the renormalization of the coupling g is given by 

y r7— n/2 

9Bare — ^T^2 9i 

and the product (16. 9p is renormalized as 

9^ = 9Bare^ Bare- 

Hence, one has the same equation as (16. 2p with solution (16. 7p but with = Z'^^'^ and 
7r = — n/272. (Recall that the anomalous dimension 72 is defined with respect to the 
renormalization constant 2'^^.) 

Furthermore, one can construct the so-called invariant charge by multiplying the product 
(16. 9p by the corresponding propagators 



i = 9T(^,9]f{D'''{%9]. (6.10) 



The invariant charge ^, being RG-invariant, obeys the RG equation without the anomalous 
dimension and plays an important role in the formulation of the renormalization group 



en 



together with the effective charge. In some cases, for instance in the MOM subtraction 
scheme, the effective and invariant charges coincide. 

The usefulness of solution (16.71) is that it allows one to sum up an infinite series of logs 
coming from the Feynman diagrams in the infrared {t —>■ —oo) or ultraviolet {t oo) 
regime and improve the usual perturbation theory expansions. This in its turn extends the 
applicability of perturbation theory and allows one to study the infrared or the ultraviolet 
asymptotics of the Green functions. 

To demonstrate the power of the RG, let us consider the invariant charge in a theory with 
a single coupling and restrict ourselves to the massless case. Let the perturbative expansion 
be 

^{^,g) = g{l + bgln^ + ...). (6.11) 

The f3 function in the one-loop approximation is given by 

Pig)=bg'. (6.12) 

Notice that the coefficient b of the logarithm in eq. (16.111) coincides with that of the (3 function. 
Alternatively the (3 function can be defined as the derivative of the invariant charge with 
respect to logarithm of momentum 

P{9)=p'^,a^,,9)\p^=,^. (6.13) 

This definition is useful in the MOM scheme where the mass is not considered as a 
coupling but as a parameter and the renormalization constants depend on it. We will come 
back to the discussion of this question below when considering different definitions of the 
mass. 

According to eq. (l6.7p (with vanishing anomalous dimension) the RG- improved expression 
for the invariant charge corresponding to the perturbative expression (16. lip is: 

^RG{^,g) = ^PT{i,g{^,g)) = g{^,g), (6.14) 

where we have put in eq. (16.71) = fi^ and then replaced t hj t = Inp^/yU^. The effective 
coupling is a solution of the characteristic equation 

-git,g) = bf, giO,g) = g, t = \n^. (6.15) 
at /i^ 

The solution of this equation is 

S(M) = ^. (6.16) 

Being expanded over t, the geometrical progression (16.161) reproduces the expansion (I6.1ip : 
however, it sums the infinite series of terms of the form (/"t". This is called the leading log 
approximation (LLA) in QFT. To get the correction to the LLA, one has to consider the 
next term in the expansion of the /3 function. Then one can sum up the next series of terms 
of the form g^f^'^ which is called the next to leading log approximation (NLLA), etc. This 
procedure allows one to describe the leading asymptotics of the Green functions for t ±oo. 
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Consider now the Green function with non-zero anomalous dimension. Let its perturba- 
tive expansion be 

T{^,g) = l + cg\n^ + ... (6.17) 

Then in the one-loop approximation the anomalous dimension is 

l{g) = eg. (6.18) 

Again the coefficient of the logarithm coincides with that of the anomalous dimension. In 
analogy with eq. fl6.13p the anomalous dimension can be defined as a derivative with respect 
to the logarithm of momentum 



7(<7)=p^^lnr(^,^?)|,.=,.. (6.19) 

Substituting fl6.18p into eg. (16. 71) . one has in the exponent 

t g , . 9 

I^i9i^^9)dt = lj^^dg = l^dg = l\n'-. 
g g ^ 

This gives for the Green function the improved expression 

l + ct + ... (6.20) 
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Thus, one again reproduces the perturbative expansion, but expression (16.201) again contains 
the whole infinite sum of the leading logs. To get the NLLA, one has to take into account 
the next term in eg. (16. 181) together with the next term of expansion of the (3 function. 
All the formulas can be easily generalized to the case of multiple couplings and masses. 



6.2 The effective coupling 

By virtue of the central role played by the effective coupling in RG formulas, consider it 
in more detail. The behaviour of the effective coupling is determined by the /3 function. 
Qualitatively, the f3 function can exhibit the behaviour shown in Figj23l We restrict ourselves 
to the region of small couplings. 

In the first case, the /3-function is positive. Hence, with increasing momentum the effec- 
tive coupling unboundedly increases. This situation is typical of most of the models of QFT 
in the one-loop approximation when i3{g) = bg^ and 6 > 0. The solution of the RG eguation 
for the effective coupling in this case has the form of a geometric progression (16.161) . It is 
characterized by the presence of a pole at high energies, called the Landau pole. We will 
consider this pole in detail later. 

In the second case, the /3-function is negative and, hence, the effective coupling decreases 
with increasing momentum. This situation appears in the one-loop approximation when 
6 < 0, which takes place in the gauge theories. Here we also have a pole but in the infrared 
region. 

In the third case, the /^-function has zero: at first, it is positive and then is negative. 
This means that for small initial values the effective coupling increases; and for large ones. 
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Figure 23: The possible form of the /5-function. The arrows show the behaviour of the 
effective couphng in the ultraviolet regime (t ^ oo) 



decreases. In both the cases, with increasing momentum it tends to the fixed value defined 
by the zero of the /^-function. This is the so-called ultraviolet stable fixed point. It appears 
in some models in higher orders of perturbation theory. 

Eventually, in the last case one also has the fixed point but now for the small initial 
coupling it decreases and for the large one it increases, i.e., with increasing momentum the 
effective coupling moves away from the fixed point, it is ultraviolet unstable. On the contrary, 
with decreasing momentum it tends to the fixed point, i.e., it is infrared stable. It appears 
in some models in lower dimensions, for instance, in statistical physics. 



6.3 Dimensional regular izat ion and the MS scheme 

Consider now the calculation of the (3 function and the anomalous dimensions in some 
particular models within the dimensional regularization and the minimal subtraction scheme. 
Note that in transition from dimension 4 to A — 26 the dimension of the coupling is changed 
and the "bare" coupling acquires the dimension [qb] = 2e. That is why the relation between 
the "bare" and renormalized coupling contains the factor (/i^)^ 

qb = ifi'rZgg. (6.21) 

Hence, even before the renormalization when Zg = 1, in order to compensate this factor the 
dimensionless coupling g should depend on fi. Differentiating (I6.2ip with respect to /i^ one 
gets 

dlogZg dg 
= eZgg + ——^Zgg + Zg 



dlogn'^ ^ ^ d log n'^ ' 
i.e., 

^ / X dq dlogZa ^ , ^ , , 



In the MS scheme the renormalization constants are given by the pole terms in 1/e 
expansion and so does the bare coupling. They can be written as 



1 + E^ = 1 + EE^- (6.23) 



n=l ^ n=\ m=n ^ 
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n=l m=n 



(6.24) 



Differentiating eg. (16.231) with respect to In/x^ and having in mind the definitions (16.51) 
and (16. 6p . one has: 

-[1 + E ^]7r(^7) = [-eg + P{9)] | E 



71=1 



n=l 



Equalizing the coefficients of equal powers of e, one finds 



7r(^) = g-^ci{g), 
dg 



g-i-Cn{g) 

dg 



d 



br{g)+f3{g)^]cn-i{g), n>2. 



(6.25) 
(6.26) 



One sees that the coefficients of higher poles c^, n > 2 are completely defined by that 
of the lowest pole ci and the (3 function. In its turn the /3-function is also defined by the 
lowest pole. To see this, consider eq. (16.241) . Differentiating it with respect to In/x^ one has 



9 + E 



a-nig) 



71 = 1 



+ [-eg + (3{g)] 



dgh ^" 



0. 



(6.27) 



Equalizing the coefficients of equal powers of e, one finds 



d 



Kg - ^^""^^^ 



{g^ - l)ai(^), 
dg 

Pig)-^an-i{g), n>2. 



(6.28) 
(6.29) 



Thus, knowing the coefficients of the lower poles one can reproduce all the higher order 
divergences. This means that they are not independent, all the information about them is 
connected in the lowest pole. In particular, substituting in (16.291) the perturbative expansion 
(I6.24P one can solve the recurrent equation and find for the highest pole term 
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i.e. in the leading order one has the geometric progression 

.2e 9 



gB = fj' 



1 - gau/e'' 



(6.30) 



(6.31) 



which reflects the fact that the effective coupling in the LLA is also given by a geometric 
progression (16.161) . 

The pole equations are easily generalized for the multiple couplings case, the higher poles 
are also expressed through the lower ones though the solutions of the RG equations are more 
complicated. 

Consider now some particular models and calculate the corresponding /^-functions and 
the anomalous dimensions. 
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The 4>^ theory 



The renormahzation constants in the MS scheme up to two loops are given by eqs. 
fi4.lOy4.i4y4.l8p . {g = A/167r2) 

Z. - l + + (6.32) 

V = 1 + ^' (6-33) 
Zn = l + —g + g\A,-—). (6.34) 

Notice that the higher pole coefficient 022 = 9/4 in the last expression is the square of the 
lowest pole one an = 3/2 in accordance with eq. fl6.30p . 
Applying now eqs. (l6.25p and (16.280 one gets 

74(^7) = ^g-^g', (6.35) 
i2{g) = ^g\ (6.36) 

P{g) = ^(^, + 272) = ^/-^/. (6.37) 

One can see from eq. (16.370 that the first coefficient of the /3-function is 3/2, i.e., the 0^ theory 
belongs to the type of theories shown in Figj23k). In the leading log approximation (LLA) 
one has a Landau pole behaviour. In the two-loop approximation (NLLA) the /3-function 
gets a non-trivial zero and the effective coupling possesses an UV fixed point like the one 
shown in Fig j23ll . However, this fixed point is unstable with respect to higher orders and is 
not reliable. Here we encounter the problem of divergence of perturbation series in quantum 
field theory, they are the so-called asymptotic series which have a zero radius of convergence. 

QED 

In QED in the one-loop approximation the renormahzation constants in the Feynman 
gauge are given by eq. (15.20 . Due to the Ward identities the renormahzation of the coupling 
is defined by the photon wave function renormahzation constant and is gauge invariant. 
Equation (15. 2p allows one to determine the anomalous dimensions and the /3-function 

7i(a) = —a, (6.38) 
72(a) = a, (6.39) 

73(a) = ^a, (6.40) 
7^(a) = -4a, (6.41) 
Paicy) = (6.42) 
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where we use the notation a = e^/16 

Thus, in QED in the one-loop approximation the effective coupling behaves the same 
way a in the 0^ theory and has a Landau pole in the LLA. In this theory, the next term of 
expansion of the /3-function is also calculated. It has the same sign. 



QCD 

In QCD the calculation of the (3 function can be based on various vertices. The result 
should be the same due to the gauge invariance. To simplify the calculations, we choose the 
ghost-ghost-vector vertex. The renormalization constants in the one-loop approximation in 
the Feynman gauge are given by (15. 8p and lead to the following anomalous dimensions and 
the /3-function: 

C 

7i(tt) = — (6-43) 
C2 

72(a) = -ya, (6.44) 
73(a) = -{^C2-'^nf)a, (6.45) 
/3„(a) = a(27l + 272 + 73) = -(yC-^r^;)a^ (6.46) 

where like in QED we take a = g'^/lGir'^, the Casimir operator C in the case of SU(3) groups 
is equal to 3, and Uf is the number of quark flavours. 

One can see from eq. fl6.46p that if the number of flavours is less than ^6*2 = Y' 
/3-function is negative and the effective coupling decreases and tends to zero with increasing 
momentum. This type of behaviour of the effective coupling is called the asymptotic freedom. 
It takes place only in gauge theories. 



6.4 Aqcd 

The solution of the characteristic equation for the effective coupling, which is a differential 
equation of the first order, depends on initial conditions. Therefore, the solution 06.161) 
depends on the choice of the initial point and the value of the coupling at this point. However, 
this choice is not unique and one can choose another initial point and another value of the 
coupling and still get the same solution, as it is shown in FigJJH 

9(t,g) 




Figure 24: Different parametrizations of the effective coupling. Each curve is characterized 
by a single parameter A 



In fact, every curve is not characterized by two numbers (the initial point and the cou- 
pling), but by one number and the transition from one curve to another is defined by the 



change of this number. To see this, consider the one-loop expression for the effective couphng 
in a gauge theory and rewrite it in equivalent form 



where we have introduced the notation 



A2 = ^2g /3oa^_ (g48) 

This quantity is called Aqcd in quantum chromodynamics and can be introduced in any 
model. The numerical value of A is defined from experiment. 

Equation fl6.48p can be generalized to any number of loops. For this purpose, let us 
rewrite the RG equation for the effective coupling in the Gell-Mann - Low form. One has 

ta^=r^. (6.49) 



Combining the lower limit with In /i^ one gets 



A2 J (3,{gy 

where 



'.>S-r^. (6.50) 



A^ = ^?exp [ ) , (6.51) 



which is the generalization of eq. fl6.48p for an arbitrary number of loops. 

The quantity A, introduced this way, is //-independent but depends on the renormaliza- 
tion scheme due to the scheme dependence of the /^-function. However, the scheme depen- 
dence of A is given exactly (!) in one-loop order. Indeed, since A does not depend on /i, let 
us choose IX in such a way that ^ 0. Then for the /3-function one can use the perturbative 
expansion 

Pa{a) = Poa'^ + Pia^ + ... 



f da 1 , ^, , 



or 

da 1 

P{a) Poa 

In this limit the ratio of two parameters A belonging to two different schemes is 

A? 



In- 

^2 



A2 
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■ 1 1 ■ 
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Po 



[ci - C2] , (6.52) 



where the coefficients ci and C2 are calculated in the one-loop order. They can be found 
from perturbative expansion of any physical quantity in two different schemes 

R = gi{l + cigi + ...) 

= 5-2(1 + 026-2 + •••)• 



Since A does not depend on g, one can take any value of and eg. (16.521) is always valid. 
The difference ci — C2 does not depend on a particular choice of R (though each of them 
depends) and is universal. 

It should be noted that the quantities like the invariant or effective coupling, the (3- 
function, etc. are not directly observable. Therefore, their dependence on the subtraction 
scheme does not contradict the independence of predictions of the method of calculations. 
We perform the perturbative expansion over the coupling which is scheme dependent, but 
the coefficients are also scheme dependent. As a result, within the given accuracy defined 
by the order of perturbation theory the answer is universal. 

In the minimal subtraction schemes when the renormalizations depend only on dimension- 
less couplings, the one-loop renormalization constants and hence the anomalous dimensions 
and the /^-function are the same in all schemes; the difference starts from two loops. The 
exception is the /3-function in a theory with a single coupling like QED, QCD or the 0^ the- 
ory, where the difference starts from three loops. Indeed, if one has two subtraction schemes 
Ml and M2 so that the couplings in two schemes are related by 

92 = q{9i) =9i + cgl + 0{gl), 
then the /^-functions (3i{gi) and P2{92) are connected by the relation 

dgi 

and their perturbative expansions are 

p2{g2) = Pog2 + Pigl + P2gl + -- 

so that the first two terms of the /3-function are universal. 

As for the further terms of expansion, they depend on the renormalization scheme and 
one can use this dependence as discretion, for instance, one can put all of them equal to 
zero. Then we would have an exact /5-function. However, one should have in mind that it 
is not valuable by itself but rather in the aggregate with the PT expansion for the Green 
functions for which we construct the solution of the RG equation. This expansion in our 
"exact" scheme is unknown. 

6.5 The running masses 

In the minimal subtraction scheme the renormalization of the mass is performed the same 
way as the renormalization of the couplings, i.e., the mass is treated as an additional coupling 
and is renormalized multiplicatively, namely, 

where the mass renormalization constant is independent of the mass parameters and 
depends only on dimensionless couplings. Then, in full analogy with the effective coupling 
one can introduce the effective or the "running" mass 

^Th{t,g) = rfi'ymig), m{0,g)=mo. (6.53) 
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Solving this equation together with the equation for the effective couphng f l6.8p one has 



m(t, g) = moG o 



/ 



lm{g{t,g))dt 




(6.54) 



In the one-loop order 



7m(a) = ca 



and the solution is 




) 



c/b 



m{t) = mo 



This is the running mass! 

The natural question arises: what is the physical mass measured in experiment and how 
is it related to the running mass and at what scale? 

To answer this question, consider why the mass is running. This is due to the radiative 
corrections. If one considers the value of momentum which is bigger than the mass, i.e. > 
m^, then the particles are created, they are running inside the loops and give the contribution 
to the running. On the contrary, if p"^ < m^, particles are not created, they "decouple" and 
do not contribute to the running. In the MOM scheme this takes place automatically because 
for the momentum smaller than the mass the diagram simply disappears. In the minimal 
scheme, on the contrary, this does not happen. Hence, it is quite natural in this case to stop 
the running at the value of = and to identify the physical mass with the running mass 
at the scale of the mass, i.e 



However, this is true only up to finite corrections. Let us come back to the definition 
of the mass term in the Lagrangian. It is chosen in such a way that the propagator of a 
particle, which is the inverse to the quadratic form, has the pole at = m^. Therefore, a 
more appropriate definition of the physical mass is the position of the pole of the propagator 
with allowance for the radiative corrections, .i.e.. 



This definition of a mass does not depend on a scale and it is also scheme independent 
and may have physical meaning. The pole mass can be expressed through the running mass 
at the scale of a mass with finite and calculable corrections. 

Consider as an example the quark mass in QCD. The quark propagator is graphically 
presented in FigJ2Sl 




physical mass = pole mass 



+ 




+ 



Figure 25: The quark propagator 



The corresponding expression is 
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The pole mass is now defined as a root of the equation 

p(l + _ rn(l - B{p^)) = 

which gives in the lowest order 

1 - B(m^ 



(6.55) 
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1 + A(m2) 

To calculate the functions A and i?, consider the one-loop diagram shown in Figl26l 




p-k 



Figure 26: The quark propagator in one loop in QCD 



The corresponding expression is 

gl „ f dk 'j^'{p - k + m)Y g^" 
~~(2^ [{p-ky-m^] "fc2" 

and was calculated earlier. The result has the form (13.191) 
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After subtraction of divergences in the MS'-scheme one has 
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A^^{p\m^) 
B'^{p',m^) 
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p^x{l — x) — m? 



Substituting p^ = w?, one finds 



A^^(m^ m^) = 2 + In ^, ^^^^(m^, m^) = -6 - 4 In 



(6.56) 
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Thus, for the radiative correction to the pole mass we have 



l + ;!i£f:(4 + 31n4) 



(6.62) 



Substituting Cp = 4/3 and ^'^ = one obtains the desired relation between the pole mass 
and the running mass at the mass scale 



TTipoie = m{m) 



4q;3 

3 TT . 



(6.63) 
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7 Lecture VII: Zero Charge and Asymptotic Freedom 



Since the behaviour of the effective couphng has so essential consequences we consider two 
typical examples which are realized in quantum field theory in the one-loop approximation 
and presumably take place in a full theory. Usually, one speaks about the zero charge 
behaviour or the asymptotic freedom. We explain below what it means. 



7.1 The zero charge 

The notion of the zero charge appeared in QED in the leading log approximation. This is 
what takes place within the renormalization group method in the one-loop approximation. 
If one writes down the expression for the renormalized coupling as a function of the "bare" 
coupling, i.e. inverts eq. fl6.3ip . one gets 

^ 9b ^ 9b 

^~ l + (3o9B/e ~ l + /?o^BlogA2' ^ ' > 

where the first coefficient of the /^-function /?o > 0. Then, removing the regularization, i.e., 
for £ — » or A — s> oo, the renormalized coupling tends to zero independently of the value 
of the "bare" coupling. This is what is called the zero charge. For the effective coupling 
considered above the zero charge corresponds to the behaviour shown on the left panel of 
Figj27] which is characterized by the Landau pole at high energies. 




Figure 27: The behaviour of the effective coupling: the zero charge (left) and the asymptotic 
freedom (right) 



The zero charge behaviour is typical of QED, the 0^ theory for positive quartic coupling 
and also the Yukawa type interactions, i.e., in those theories where the /^-function is positive. 
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It is obvious that in the vicinity of the pole the perturbation theory does not work and, 
hence, the one-loop formula is not applicable. However, for small momenta transfer the 
one-loop approximation is reliable. For instance, in QED the effective expansion parameter 
is e^/167r^ = a/Air ^ l/137/47r ^ 5.8 ■ 10~^ and the next loop corrections (which have the 
same sign) do not play any essential role. The behaviour of the effective coupling in QED 
in the region up to 100 GeV has got the experimental confirmation in measuring the fine 
structure constant at the LEP accelerator. At the scale equal to the mass of the Z-boson Mz 
the fine structure constant is not 1/137 but a{Mz) ~ 1/128, which is in a good agreement 
with the one-loop formula. 

The large momenta transfer in this case are limited by the pole provided the pole does 
not disappear in a full theory. It is still unclear how higher orders of perturbation theory 
influence this behaviour since the perturbation series is divergent and it is impossible to 
make definite conclusions without additional nonperturbative information. 

The presence of the Landau pole indicates the presence of unphysical ghost states. To 
see this, consider the photon propagator in QED which due to the Ward identities coincides 
with the invariant charge and in the leading log approximation has the form of a geometric 
progression 



where Q is the electric charge of a particle (in the units of electron charge) running round 
the loop. 

This expression has a pole in the Euclidean region at = — m^ea:p( ^^Q2^ ). Substituting 

m = me = 0.5 MeV, ao ~ 1/137 and EQ^ = [(4/9 + 1/9)3 + 1)3] = 8, one gets ~ 
— (5 ■ 10'^^)^ GeV^. That is the pole is very far off, even beyond the Planck scale, and at low 
energies one can ignore it. However, the presence of the pole indicates the presence of a new 
asymptotic state and the residue at the pole defines the norm of this state. In the case of 
the Landau pole the residue is negative, i.e., the new state is a ghost, it has the wrong sign 
of the kinetic term in the Lagrangian. This fact, in its turn, leads to negative probabilities, 
which indicates internal inconsistency of the theory. 

Usually, it is assumed that there are two ways out of this trouble: either the higher order 
corrections improve the behaviour of the theory at high momenta so that the Landau pole 
disappears, or that the zero charge theory is contradictory by itself, but at high energies 
it is part of a more general theory where the behaviour of the coupling is improved. The 
example of such a behaviour is given by the Grand Unified Theories where QED is one of 
the branches of a non-Abelian gauge theory with the asymptotically free behaviour. In both 
the cases the theory at high energies is modified. At the same time, the zero charge theory 
is infrared free, i.e. for small momenta transfer the coupling goes to zero. 

7.2 The asymptotic freedom 

The name asymptotic freedom originates from the non-Abelian gauge theories where it was 
found that the sign of the first coefficient of the /3-function is negative. The effective coupling 
in this case behaves as is shown in the right panel of FigJ27] and tends to zero at high 
momenta transfer. This means that quarks in QCD are quasi- free particles, i.e., practically 
do not interact. This way one explains the success of the so-called parton model of the 
strong interactions at high energies, according to which the proton behaves as a set of free 




(7.2) 
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partons, and at high energies the interaction takes place with the individual partons and 
their interaction does not play any role. 

The behaviour of the effective coupling in QCD at high energies was tested at various 
accelerators and in various experiments and the validity of the renormalization group formula 
was confirmed. The accuracy of modern measurements assumes the inclusion of the next 
terms of perturbative expansion. In QCD in the MS scheme the four terms of the /3-function 
are known. Below we present the two-loop expression 

= -^[11 - InfK - (i^[102 - ^nf]al + O(a^). (7.3) 

As one can see, if the number of quarks in not too big, both the coefficients of the /3-function 
are negative. All the experimental data fit a single curve for the effective coupling with the 
parameter Kqcd ^ 200 MeV (see Figj28D 




Figure 28: The variation of the effective coupling of the strong interactions as with energy 

In four-dimensional space the asymptotic freedom occurs only in non-Abelian gauge 
theories. But in the case when one has several interactions, like in the Standard Model, 
the non-Abelian coupling may draw other couplings into the asymptotically free region. 
Consider, for instance, the behaviour of the Yukawa couplings in the SM. For simplicity, let 
us take a single Yukawa coupling for the t-quark and a single gauge coupling. Then in the 
one-loop approximation the equations for the effective couplings look like 



where the coefficients b, a and c are always positive and for the SM are equal to 7, 9/2 and 
8, respectively. The solutions to these equations are 
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(7.5) 



E{t) = ig/go f\ Fit) = f E{t')dt'. 
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In the case of a single Yukawa coupling it can be written in an explicit form 
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1/0 1 



30 ' 



(7.6) 



l + ^^[(^)=^-i-ll' 

go c-b go ^ J 

Graphically, it can be presented in a phase diagram shown in Figj29l For the initial condition 




9 

Figure 29: The behaviour of the Yukawa and gauge couplings for various initial conditions 

such that yo > {c — b)/agQ the Yukawa coupling increases with momenta and has the Landau 
pole, while for ?/o < {c — b) / a go it demonstrates the asymptotically free behaviour. In a 
similar way in the Grand Unified Theories one can reach the asymptotic freedom for all the 
couplings. 

The back side of the asymptotic freedom at high energies is the presence of a pole at low 
energies or the infrared pole. In this region, we also go beyond the validity of perturbation 
theory since the coupling increases. To find the true behaviour of the coupling one has to 
attract independent nonperturbative information. However, in QCD the region near the 
infrared pole p ~ Aqcd is in the phase of hadronization, i.e., in this region the quark-gluon 
description is no more adequate. Therefore, the behaviour of the effective coupling in this 
region is not described by perturbative QCD. 



7.3 The screening and anti-screening of the charge 

The variation of the coupling with momenta transfer or with the scale, which is the char- 
acteristic feature of quantum field theory, has its analog in a classical theory. This analogy 
allows one to understand the qualitative reason for the variation of the coupling. 

Indeed, let us consider the electromagnetic phenomena. Consider the dielectric medium 
and put the test electric charge in it. The medium will be polarized. The electric dipoles 
present in the medium will be rearranged in such a way as to screen the charge (see Fig j30|) . 
This is a consequence of the Coulomb law: the opposite charges are attracted and the same 
charges are repulsed. This is the essence of the electric screening phenomena. 
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Figure 30: The electric screening and magnetic anti-screening 



The opposite situation occurs in magnetic medium. According to the Bio-Savart law, the 
electric currents of the same direction are attracted and the opposite direction are repulsed 
(see Fig j30l) . This leads to the anti-screening in magnetic medium. 

In quantum field theory the role of the medium is played by the vacuum. The vacuum is 
polarized in the presence of created virtual pairs. The matter particles as well as transversely 
polarized quanta of the gauge fields act like the electric dipoles in the dielectric and cause the 
screening of the charge. At the same time, the longitudinal quanta of the gauge fields behave 
like currents and cause the anti-screening. These two effects are in competition (see eq. (13.311) 
above) and, for instance, in QCD with a small number of quarks the effect of anti-screening 
prevails. 

Thus, the couplings become the functions of the distance or momentum transfer described 
by the renormalization group equations. 



8 Lecture VIII: Anomalies 



The gauge invariance leads to numerous relations between various operators and their vac- 
uum averages, i.e., the Green functions. We have already come across such relations called 
the Ward or the Slavnov- Taylor identities. They are the consequences of the gauge symme- 
try of the classical theory. In case when one has divergences in a theory and is bound to use 
some regularization, the validity of these identities depends on invariance of the regulariza- 
tion. However, one can always perform the subtraction of divergences in such a way that 
the finite parts obey these relations. 

The exception from this rule is the so-called anomalies. By anomalies one usually means 
the violation in quantum theory of some relation, for instance, the conservation of the current 
or the Ward identity following from the symmetry properties of a classical theory. The well- 
known examples of quantum anomalies is the anomaly of the trace of the energy-momentum 
tensor or the axial anomaly. The characteristic feature of the anomaly is the impossibility 
of its removing by the redefinition of any quantities or parameters. 

8.1 The axial anomaly 

Consider quantum electrodynamics. Let us define the vector and the axial vector currents 

j^ = tPY^, jJ = V;7W'. (8.1) 

In classical theory the equations of motion lead to the conservation or partial conservation 
of the current 

^mJ/. = 0, d,jl = 2^mJ^ (8.2) 

where = ip'-f^ip- 

On the other hand, as a consequence of the gauge invariance, the vector and the axial 
vertices obey the Ward identities 

{p-p'n,{p,p') = S-\p) ~ S-\p'), (8.3) 
ip-p'nlip,/) = S-\p)j' + j'S~\p') + 2mT\p,p'), (8.4) 

where F^, and F^ are the vector, axial and pseudoscalar vertices, respectively, and 5* is 
the fermion propagator. 

If one looks how the identities fl8.3|8.4l) are fulfilled in perturbation theory, one first of all 
has to introduce some regularization due to the presence of the ultraviolet divergences. If the 
regularization is gauge invariant, then the vector Ward identity is satisfied in any order of 
FT. For the axial identity there are two types of diagrams: in the first one the axial current 
is in the outgoing fermion line, and in the second one the axial current is in the internal loop 
(see FigJSTl). For the first type of a diagram the identity (18. 4p is satisfied, and for the second 
type there exists one famous triangle diagram (see Fig j32|) where it is violated due to the 
ultraviolet divergence of the integral. 



Figure 31: The diagrams with the axial current in external and internal fermion lines 




Figure 32: The anomalous triangle diagram for the axial current 



Indeed, the corresponding integral in momentum space looks like 
p 

{-){-ief f -^Tr 



rw V 



(277) 



^ ^ k^^ {p + ky^ {q + kf 



.5) 



q-p 

and is formally divergent requiring the regularization. 

To preserve the conservation of the gauge invariance, it is useful to introduce the dimen- 
sional regularization; however, here we for the first time face a problem since the 7^ matrix 
has no natural and consistent continuation to non-integer dimension. Two properties of the 
7^ matrix, namely, the anticommutation with all 7^^, /x = 0, 1,2,3 and the property of the 
trace Tr(7^7'^7''7''7'^) = — 4ie^'^''°" are in contradiction if the dimension is noninteger. To 
calculate the axial anomaly, we use the following trick: we use the formula for the trace but 
reject the property of anticommutativity of 7^. This allows one to perform al the calculations 
in a consistent and unambiguous way. 

The divergence of the axial current can be obtained by multiplication of (18.51) by iq^ 
which gives 



d'^k Tr q'y^kY{k + p)'J^{k + q) 



(277) 



F(A; + p)2(fc + g)2 



1.6) 



Using the cyclic property of the trace we move q to the right and write it as g = (g + A;) — k. 
Then the first term multiplied by A; + g gives {k + g)^ and cancels with the denominator. As 
a result, one gets the integral 

d^k TrUj^kYik + pW 



(277)4 k^{k + py 

which depends only on p and after the integration turns to zero due to the antisymmetry of 
the trace with the 7^ matrix. 
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In the second term we will drag k to the left until it is multiplied by k giving k'^. As a 
result, at each step we always get the trace of four 7-matrices with 7^ for which we have the 
formula with the e-tensor. We obtain in the numerator 

_4 •^ai./3A^a^^ + + q)^ - q^] + Sie'^'^^Pk^'ik + pfq^k^ - Aie'"'^Pk''q%k + pf - p^] 

_^-^uaXppa^Pf^2 ^ sie^'^^^k^'p'^ q^k^ 

Despite the fact that the integral is formally divergent, using a dimensional regularization 
and collecting all terms together we finally get the finite answer equal to 

p2 2 



r-^rf 



= T^e^^^'^M^.^p^A = -—e^^'^F^^F^^. (8.8) 



One has to add to this expression the same diagram but with the replacement p ^ 
q — PiV ^ A and take the sum, but the answer is already invariant with respect to this 
replacement. Multiplying (18.71) by Ay{j>)Ax{{q — p) and transforming to the coordinate 
representation, one gets 

As a result one has the following modification of equations for the divergence of the axial 
current and the axial vertex 

d,3l = 2tmf + £F^,F,.e^'^^^ (8.9) 

{p-p'nUP^P) = S'\ph' + l'S-\p') + 2mT\p,p') - i—F{p,p'), (8.10) 

where F{p,p') is the vertex with insertion of the operator FF. The appearance of the r.h.s 
in these equations is called anomaly known as the Adler-Bell-Jackiw or triangle anomaly. 

The most essential here is not the violation of the Ward identity but the fact that sub- 
tracting the anomaly and restoring the "normal" Ward identity for the axial vertex we 
violate the conservation of the vector current. In other words, it is impossible to satisfy the 
conservation of axial and vector currents simultaneously. 

Notice that the violation of the conservation of the axial current preserving the con- 
servation of the vector current (18.91) can be obtained by accurately calculating the matrix 
element for the divergence of the axial current in x-space splitting the arguments of the field 
operators. Consider the vacuum average of the divergence of the axial current, and to avoid 
the singularity for the product of two operators at coinciding points, split the arguments. 
Then to preserve the gauge invariance, we have to insert between the operators the exponent 
of the Wilson line. The axial current then takes the form 

x+e/2 

jj(a;) = lim{Vi(a: + e/2)7Vexp[-ie J dz" A^{z)]ij{x - e/2)}, (8.11) 

x-e/2 

and for the divergence we get 

/'x-\-e /2 

d.f(x) = lim{a.V'(a; + 5/2)7 Vexp He / dz" AJz)U{x - 6/2) 

^ Jx-e/2 

px-\-e /2 

+ ip{x + e/2)Yl^exp[-ie dz"" A^{z)]dMx - e/2) (8.12) 

Jx-e/2 

/*x-\-e /2 

+ '4!{x + e/2)-f^-f^[-iee''d.A^{x)]exp\-ie dz''A^{z)U{x-e/2)}. 

Jx-e/2 



Using the equations of motion 

l^d^ip = —ieAip, d^ip'y'^ = ieipA 
and keeping the terms of the order of e we find 

d^il^x) = hm{(9^^(x + e/2)[-ieA{x + e/2) - ieA{x - e/2) 

-tee''rd^A,{x)]^''^lj{x - e/2)} 
= \im{^{x + e/2)[-tee''Y{df,A, - d,A^)]'y''^{x - e/2)} (8.13) 

Now we have to calculate the vacuum average over the fermion vacuum (the photon field is 
assumed to be external) which means that we have to permute the fermion operators. The 
permutation function of the fermion operators is singular and this is the reason for appear- 
ance of a nonzero term similarly to the appearance of triangle anomaly due to divergency of 
the integral. Indeed, calculating the propagator of the fermion in external field and keeping 
the terms linear in the photon field, we get 

The propagator fl8.14p is singular as y ^ z; however, the first term does not give a contri- 
bution to the divergence, while the second one leads to 

{ip{x + e/2)Yl^i){x - e/2)) = 



Tr[^;r^,HeA{p))^.^^'] 



(27r)4(27r)4 '{k + p)^' '"^''P 

d'k d'p ^,^,^_,,^ 4ee^'^P-{k + p),A,{p)k^ ^ 



(27r)4(27r)4 {k + p)^k^ 

To find the limit as £ — >^ 0, one can expand the integrand for large k, which gives 

(V5(x + e/2)7VV^(x-e/2)) = 4ee^-P- / -^e^>.A,(p) /^e'^'^^^ 



(27r)4 ^^^^y (27r)4 k^ 

-4ee^-^-.9.A,(x)^ = -ee'^''^^F.,{x)^„ (8.16) 



Substituting this expression into f l8.13p we find 



d,jl = lim{-e6'^^''-2F,,(x)^(-^ee-F^.)} = -^e^''^'^ F,,F,„ (8.17) 

that coincides with (18. 9p . 

The axial anomaly has one very important property: the obtained formulas (18. 9p and 
(18.1 op are exact in all orders of perturbation theory, i.e., have no radiative corrections. More 
rigorous statement is: there exists such a renormalization scheme (and it was constructed 
explicitly) that the radiative corrections to the axial anomaly are absent. This statement is 
the subject of the Adler-Bardeen theorem. Graphically, this means the cancellation of the 
contributions of the diagrams shown in Figj33| which was checked by explicit calculation. 

The Adler-Bardeen theorem is valid also in non-Abelian theories. It has important con- 
sequences: if the anomaly is compensated in the lowest order, it will not appear further. 




Figure 33: Cancellation of radiative corrections to the axial anomaly 
8.2 Consequences of the axial anomaly 

Let us ask the question what are the consequences of the axial anomaly? Here one has to 
distinguish two cases: when the operator of the axial current is an external operator with 
respect to the Lagrangian and when it is present in the interaction Lagrangian. 

In the first case, the presence of anomaly does not lead to any troubles and even may 
be useful. Thus, for instance, in the current algebra which describes the low energy hadron 
interactions, the axial anomaly is responsible for the neutral pion decay vr" — ^ 27 and is in 
agreement with the experiment. 

In the second case, the triangle anomaly leads to that the ultraviolet renormalizations of 
the vector vertex do not remove all divergences from the axial vertex. This has destructive 
consequences for the renormalizability of the whole theory. To see this, compare the two 
processes of the elastic scattering of leptons: i/g + e — i^e + e and z/^ + e — z/g + /i in the 
Standard Model. Graphically, in the lowest order they differ by one diagram containing the 
triangle anomaly (See Fig JM|l . 




Figure 34: The anomaly in the process of lepton scattering in the Standard Model 

As a result, after the renormalization the amplitude of i^^e-scattering has finite radiative 
corrections, while that of z/e-scattering is divergent. This led to nonrenormalizability of the 
theory and was a serious problem for the left-right nonsymmetric model with SUii^) x f/(l) 
symmetry before the introduction of the -quark. Remarkably, the -quark introduced by 
Glashow, Iliopoulos and Maiani for suppression of the neutral current changing strangeness 
leads to the compensation of the contributions of quarks and leptons to triangle anomaly 
and restores the renormalizability of the theory. 

In the Standard Model due to its left-right asymmetry the presence of the axial currents 
for quarks and leptons leads to several kinds of triangle anomalies where all three gauge 
fields may be in the vertices of the triangle. However, not all of them lead to anomalies. In 
general, the anomaly is proportional to the trace 

Tr T^iT^T"}, 



where the matrix T" is the generator of the corresponding gauge group in the representa- 
tion corresponding to the fields that run inside the triangle. The necessary condition of 
the existence of anomaly is the presence of the complex representations and the nontrivial 
anticommutator of the generators of the group. Among the simple Lie groups which satisfy 
this requirement, only the groups SU{n), S'0(4n + 2) and Eq have complex representations 
and out of them only the SU{n),n > 2 and SO{Q) groups have a symmetric invariant needed 
for the construction of the anomaly. The gauge theories built on other groups are free from 
anomalies. 

The non-vanishing anomalies corresponding to the symmetry group of the Standard 
Model S'f/c(3) X SUl{2) x ?7y(l) are presented in Fig 135] where the gauge fields adjusted 
to the groups U{1) and SU{2) are shown prior to mixing. The particles that run over the 
triangle can be either left or right quarks and leptons. Particles of different helicity give the 
opposite sign contribution to the axial anomaly. 




U(1) U(1) SU(2) SU(2) SU(3) SU(3) grav grav 



Figure 35: The triangle anomaly in the Standard Model 



In the first case, the anomaly is proportional to the trace of the cube of hypercharge 
of quarks and leptons in each generation 



TrYl — TrY^ and its absence is achieved by the cancellation of the contributions 



TrY^ 



= 3 




i-f- 

^3^ 


4)^- 


l--f] 


T 


T 


T 


T 


T 


colour ul 


dL 


Ur 


dR 



:-2f = 

T 

eR. 



0. 



In further diagrams the anomaly is proportional to, respectively. 
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(8.19) 



This way the anomaly is miraculously canceled in all the cases and does not break the 
renormalizability of the SM. 



8.3 The conformal anomaly 

Another example of quantum anomaly is the conformal anomaly or the anomaly of the trace 
of the energy-momentum tensor. The requirement of conformal (scale) invariance means the 
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invariance of the action with respect to the transformation 



x^, x^e"'", 0(xe"'") -> e'^'"0(x), (8.20) 

where A is the dimension of a field. This condition is fulfilled in the classical Lagrangian 
if it has no dimensional parameters. In this case, according to the Noether theorem, there 
exists a conserved current called the dilatation current = Q^^x„, so that 



where is the symmetric energy-momentum tensor. 

The easiest way to see it is to define the energy-momentum tensor as a variation of the 
action of the matter fields with respect to the space-time metric in the external gravitational 
filed 

^ 2^ / d\ C{x). (8.21) 
The scale transformation can be realized as a variation of the metric 

qA^) ^ e^'^gA^). (8.22) 

This means that the variation of the Lagrangian under this transformation is the trace of 
B^'^. The deviation of the trace of the energy- momentum tensor from zero indicates the 
violation of the scale (and hence conformal) invariance. 

In the quantum case, due the presence of the ultraviolet divergences the new scale ap- 
pears. This is the same phenomenon of dimensional transmutation discussed above. There- 
fore, the scale invariance of the action is violated. 

Since the coupling constant becomes scale dependent, its variation with the scale (18.201) 
takes the form 

5g = cTfi^ = a(3{g). (8.23) 
Hence, for the variation of the Lagrangian we get 

5 C = a^mg}), (8.24) 

dgi 

i.e., 

d.D'^ = e;= ^^ma})- (8.25) 

This relation is known as the trace anomaly of the energy-momentum tensor. 

Similarly to the axial anomaly, relation fl8.25p can be checked by perturbation theory. 
However, in this case the result is defined by the full /3-function calculated in all orders of 
PT. 
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9 Lecture IX: Infrared Divergences 

One more problem that we encounter on the way of calculating the finite expressions for the 
probabilities of physical processes is the presence of the so-called infrared divergences. They 
appear when calculating the matrix elements of the scattering matrix on shell, i.e., when the 
squares of external momenta are equal to the corresponding masses squared and the theory 
contains massless particles like photons or gluons. The infrared divergences can be of two 
types: the divergences for small values of momenta (the genuine infrared divergences) and 
the divergences at parallel momenta (the coUinear divergences). Contrary to the ultraviolet 
divergences, the infrared divergences have a clear physical meaning: a massless particle 
with a very small momentum can not be registered and with momentum parallel to another 
particle cannot be distinguished. For this reason in the theories with massless particles one 
has to define the physical process to be evaluated in a proper way. 

9.1 The double logarithmic asymptotics 

For illustration consider the process of creation of a muon pair in the e~^e~ annihilation. The 
leading diagrams for this process are shown in FiglHBl 



Figure 36: The diagrams contributing to the process e'^e~ in QED: a) the leading 

order, b)- d) the virtual corrections of the order of a, e) the real corrections of the order of 
a. 

The first diagram is the tree amplitude, it gives the contribution in the leading order. The 
radiative corrections due to emission of virtual photons (Figl36]b)) are the corrections to the 
vertex function considered above (see (13.231) ). It is easy to see that if one puts in this formula 
all fermion momenta on mass shell, i.e. = {p — qY = m?, then in the second integral in the 
denominator one gets [—m?x'^ + q^y{x — y)]. Performing the change of variables y ^ yx so 
that all the integrations are performed within the limits [0,1], we get [— m^x^ + g^x^y(l — y)], 
and the integral (with account of the Jacobian = x) is logarithmically divergent as x — ^ 0. 
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The appeared divergence has the infrared nature. Like the ultraviolet one it can be 
regularized, for instance, by introducing the nonzero photon mass or cutting the integral 
over momenta at the lower limit, or with the help of dimensional regularization but it cannot 
be removed by any renormalization. 

Let us calculate this diagram on mass shell introducing the nonzero photon mass rriph 
into the virtual photon line. This will not break the gauge invariance since, as it will be 
clear later, after the cancellation of the IR divergences one can put the mass of a photon 
equal to zero. 

Let us go back to eg. (13.231) . remove the UV divergence by the minimal subtraction and 
go to the mass shell for the fermion fields taking into account that the external fermion 
operators obey the Dirac equation [p — m)u{p) = and u{p — q){p — q — ttl) = 0. Then after 
some exercise we obtain for the vertex function the following expression: 
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where the form-factors Fi{q^) have the form 
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The form factor F2 is IR convergent and does not need any regularization. Substituting 
rriph = 0, we get 
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For = it can be easily calculated and equals 

F2{q^ = 0) = 
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(9.5) 



which is nothing else but the first correction to the g-factor, which is called the anomalous 
magnetic moment of electron (muon). 

As for the form factor Fi, it is IR divergent. We calculate its divergent part in the limit 
0. It comes only from the second integral in fl9.2p . To simplify the integration, we 
notice that the divergence is defined by the region of the parameter x ~ 0. Therefore, we 
put X = everywhere in the numerator and in the coefficient of uiph in the denominator. 
Then one gets 
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The integral over x is now easily evaluated 
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The remaining integral over y is also simple. We calculate it in the limit —q^ oo. Then 
it takes the form 

The obtained double logarithmic behaviour of the form-factor is called the Sudakov double 
logarithm. It contains the infrared cutoff in the form of the photon mass. In the amplitude 
of creation of the muon pair there are two of such form factors for the electron and the 
muon vertices, respectively. The corrections to the fermion and the photon propagators 
do not contain the IR divergences. Thus, the cross-section of the process e+e~ 
is logarithmically divergent. In order to understand the reason of appearance of the IR 
divergence and to find the method of its elimination, consider the process of creation of the 
muon pair from the point of view of an observer. 



9.2 The soft photon emission 



During the process of electron-positron annihilation the muon pair is created with momenta 
that satisfy the conservation law and can be measured. However, they are registered with 
some accuracy, and momentum smaller than some value which depends on a particular de- 
tector is not registered. Therefore, if besides the muon pair the photon with momentum 
smaller than this value is created, then this process with emission of the "soft" 7-quantum 
e^e~ —* fi^^^'-f is experimentally indistinguishable from the initial process e~^e~ — > 
The diagrams corresponding to the process e~^e~ fi~^iJ;~^ are shown in Fig 1321 e). They 
contain an additional vertex and hence additional coupling, but being squared give a cor- 
rection to the main process of the order of a, exactly as the radiative corrections due to the 
virtual photon. 

Let us compare the differential cross-sections of the precess e^e~ — > ^^fi~ in the one-loop 
approximation and e^e~ ;U"'"yU~7 in the tree approximation. We have, respectively. 
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where the second cross-section is written down without derivation which we will perform 
later. As follows from eqs. (19.9119. lOp . each of these cross-sections is IR divergent, but in the 
sum the divergences cancel and one gets the finite answer. 

What is observable after all? In fact, neither the first nor the second process is observable 
separately. In a real detector with limited sensitivity one observes the process of creation 
of the muon pair plus an arbitrary number of soft photons with the total energy below the 
sensitivity threshold. In a given order of perturbation theory we have to sum the cross- 
sections of the two processes in order to get the observed cross-section 
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The latter cross-section is given by the same formula f lQ.lOp with the replacement in the 
second logarithm of the photon energy by Emm- Thus, we get 

As one can see, for the proper statement of the problem the cross-section of the observable 
process is finite and does not depend on the IR regulator. At the same time, it depends on the 
sensitivity of the detector Emin and for improved sensitivity tends to infinity. However, this 
infinity also is not physical and is the artefact of perturbation theory: when the logarithm 
becomes large we go beyond the scope of applicability of perturbation theory and it is 
necessary to perform the summation of these corrections by analogy with what happens 
with the ultraviolet logarithms within the renormalization group method. 

Thus, the IR divergences appear due to the contributions of the photons with "soft" mo- 
menta: real with the energy smaller than Emin and virtual with momenta k"^ < What 
is important is that the momenta of fermions are on mass shell, otherwise the singularities 
in the propagator do not arise. The typical diagram of higher order contains a big amount 
of real and virtual photon lines (see Fig J57|) . 



(9.12) 



Figure 37: The hard process with creation of the soft photons 

Let us try to sum up the contributions of these soft photons. Consider first the external 
fermion line with the outgoing photons (real and virtual). 




Figure 38: The emission of the soft photons from the fermion line 
It corresponds to the following expression: 
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(9.13) 



We use now the fact that the operator u{p) obeys the Dirac equation u{p){p — m) = and 
omit the momenta fcj -C p in the numerator. Then we get 

u{p)l^^{p + m)'^^^{p + m) ■ ■ ■ = M(p)2p''i7^2(p + m) ■ ■ ■ = u{p)2p^'^2p^^ (9.14) 

Hence, eq. fl9.13p takes the form 

Hp) (e^)(e , J ■ ■ ■ (e— ^ — ). (9.15) 

The next step is the summation over all the permutations of the photon lines and the 
permutations of momenta ki. (So far we have not distinguished between the real and virtual 
photons, we will do it later.) This operation is non-trivial but leads to the simple result. 
One has 

y ±^ I = ±±...±_ (916) 

permutations P^l Pi^l + h) p{kl + ^2 + " " " + A;„) pki pk^ pkn 

The same procedure can be applied to the incoming fermion line. The difference is 
that the fermion momentum has the opposite direction which leads to the replacement of 
{p+kiY to {p—kiY in the propagator, i.e., the change of the sign p —p in the denominator. 
Collecting both factors together we get the following expression for the amplitude of emission 
of soft photons from arbitrary points of the incoming and the outgoing line (Fig j39l) : 




hard 



Figure 39: The emission of soft photons from arbitrary points of the incoming and the 
outgoing lines 



p'kx pki) \p'k2 pk2) \p'kn pkn) 

Now we have to decide which photons are real and which are virtual. The virtual photon 
can be obtained by joining the two photon momenta fcj and kj, taking k^ = —kj = k, 
multiplying by the photon propagator and integrating over k. In this way for any virtual 
photon we get the expression: 

r d'^k —ifp'p\fp' p \ 

where the factor 1/2 compensates the double counting due to permutation of ki and kj. The 
obtained integral is nothing else but the vertex function in the one-loop approximation, i.e., 
the form factor Fi{q^). 
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If the number of virtual photons equals n, one gets the product of n expressions like 
(19.181) and the factor 1/n! taking into account the permutations which do not change the 
result. The full answer is obtained with the help of summation over the soft virtual photons, 
which gives 
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(9.19) 
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At the same time, if the real photon is emitted, then instead of the propagator one has 
to multiply the amplitude by the polarization operator, sum up over all polarizations and 
integrate the square of the matrix element over the photon phase space. In this case, one 
gets the following expression: 
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which is the element of the cross-section of the process e^e~ ^'^fi~'~f. The integration 
over the modulus of the three- vector k has to be performed within the limits {mph,Emin)- 
Contracting the indices one gets 
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The first and the last integrals are equal to each other. Let us consider the last one and 
choose the frame where p=0. This gives 
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As for the second integral, we proceed in the following way: first we also choose the frame 
p = 0, and then we covariantize the answer. One has 
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Covariantizing this answer and having in mind that q = p — p' ■, p^ 
= 2m? — 2m\lp' + m? one gets 
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In the limit —q^ — oo we get the desired answer 
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coinciding with fl9.1Up . 

If there are n real photons, there are n such contributions and the symmetry factor l/n\ 
taking into account the identity of the final particles. The cross-section of the process with 
emission of an arbitrary number of photons with the energy smaller than Emin hence equals 
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Combining the results for the real and virtual photons one gets the final expression for 
the observable cross-section with emission of an arbitrary number of photons with the energy 
smaller than Ernir, 
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The obtained expression is valid in all orders of perturbation theory. The exponential factor 
does not depend on the IR cutoff but on the sensitivity of the detector. It is called the 
Sudakov form factor. When Emin tends to zero, the form factor decreases and in the limit 
Emin vanishes. This is the manifestation of the statement that he amplitude of creation 
of the fermion pair without accompanying soft photons indeed vanishes: the charged particle 
inevitably emits the low frequency electromagnetic waves. This means that the cross-section 
of elastic electron scattering without inclusion of emission of bremsstrahlung quanta should 
vanish, precisely as it follows from eq. fl9.28p . 

Let us estimate the value of the Sudakov form factor for some real process. A good 
example is the cross-section of e~^e~ annihilation into hadrons which in the leading order 
in the fine structure constant is described by one diagram with Z-boson exchange in the 
s-channel. The cross-section has a maximum in the Z-boson peak where it is described by 
the Breit-Wigner resonance formula. The energy is equal to the Z-boson mass Mz and the 
energy resolution is defined by the Z-boson width Tz- Substituting the values Mz = 91.187 



GeV, Tz = 2.496 GeV, 
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0.5 MeV, a = 1/128 into the form factor ([221 we get 
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0.648. 



As one can see, the form factor, despite the smallness of the fine structure constant, consid- 
erably departs from unity and has to be taken into account when analysing the experimental 
data. 
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9.3 The cancellation of the infrared divergences 

The considered example is typical of the QED and one can make the general statement 
concerning the infrared divergences for the elements of the S'-matrix. 

The infrared divergences in radiative corrections to the cross-section of any physical pro- 
cess in QED are cancelled in every order of perturbation theory if to the cross-section of 
the elastic process one adds the inelastic cross-section of the process with emission of an 
arbitrary number of additional photons integrated over the phase space with the requirement 
that the total photon energy does not exceed some value Emin- 

This statement is also valid for the cross-sections of the processes in non-Abelian gauge 
theories like the electroweak theory and some processes in QCD, though in this case, due 
to the self-interaction of the non-Abelian gauge fields, there is no full factorization with 
the exponentiation, and the proof of this statement presents some problem. Nevertheless, 
for many processes the result has the same form. Thus, for example, the electromagnetic 
form- factor in QCD has the same Sudakov form fl9.28p but with the replacement a — >■ Cpc^s- 

Thus, one can say that the problem of obtaining the ultraviolet and the infrared finite 
radiative corrections to the cross-sections of the physical processes is solved in two steps: first, 
with the help of the renormalization procedure one gets rid of the ultraviolet divergences, 
which is under full control in renormalizable theories; second, defining the correct physical 
process including the emission of the soft quanta, the cancellation of the infrared divergences 
takes place. 

As we will see below, this is not sufficient in non-Abelian gauge theories with massless 
gauge fields. They contain additional divergences which require some ads-inn to the described 
procedure. We will consider this question in the last lecture. 
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10 Lecture X: Collinear Divergences 



10.1 The collinear divergences in massless theory 

The obtained result fl9.28p for the cross-section of creation of the muon pair in the process 
of e"'"e~-annihilation with emission of additional soft photons is typical of the theories with 
a massive fermion and massless photons. It can be generalized to non-Abelian theories 
with massless gluon, though the gluon interactions cause some problems in proving the 
cancellation of the IR divergences. Note, however, that eq. fl9.28p contains the logarithmic 
singularity with respect to the fermion mass, and if the latter tends to zero, one has the 
new divergence. This would not cause any problem since all the fermions are massive but 
the masses of the electron and the light quarks are so small compared to the characteristic 
energies of the scattering process that with good precision it is reasonable to neglect them. As 
for the QCD, considering the processes with gluons in initial states due to the self-interaction 
of the gluons we face this problem for the gluon amplitudes. 

Let us analyse what is the reason for the appearance of the new divergence after the IR 
divergence at small photon momenta if regularized by introducing the photon mass. Consider 
for this purpose eq. (19.171) for the contribution of the real or virtual photons. The difference 
is that in one case the integration goes over the four-momentum of the virtual photon; and 
in the other case, over the three-momentum of the real photon, but what is essential that 
for the massless electron its propagator takes the form 



2pA; 2{p^k^-pk) 2{\p\\k\-\p\\k\cose) 2|p1|A^|(l - cos^) ' 

where 9 is the angle between the electron and photon momenta. (In the case of a virtual 
photon we use the fact that the contribution to the singularity comes from the region of 
photon momentum close to the mass shell.) 

Thus, the divergence appearing in the massless case comes from the integration over 
the angles and not over the modulus, as in the case of the IR divergence, and is related 
to the collinearity of momenta of two particles. For this reason it is called the collinear 
divergence. To get rid of these divergences, one can introduce the angular sensitivity of the 
detector analogously to the IR divergence. This would reflect the fact that two massless 
particles having almost parallel momenta are not distinguishable from a single particle with 
the same total momentum. Hence, the observed cross-section should include besides the 
main process the process of emission of the soft photons and the process of emission of the 
collinear photons with the kinematically allowed absolute values of momenta. 

However, in real life the quarks and leptons are massive though their masses are very 
small; therefore, the problem of collinear divergences occurs for the processes with the gluon 
fields. Since the gluons are not free particles but exist inside hadrons, any process with the 
gluons has a similar process with quarks and it is reasonable to consider them together. For 
this reason, one usually speaks about the inclusive processes where besides the particles of 
the main process one includes the creation of all kinematically allowed particles, in particular 
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the gluons. In this case, we do not impose any restriction on the gluon energy, we do not 
introduce any detector sensitivity to the energy or the angle, but sum over all the possibilities. 
It happens, however, that this is not sufficient to get the finite answer. It is necessary to 
take into account the possibility of existence of collinear gluons in the initial state, and only 
after this one can get the finite answer for the cross-section of the observable process. 

The multiloop analysis in this case is much more complicated and is the subject of the 
Kinoshita-Lee-Nauenberg theorem which states: 

The infrared and collinear divergences in a massless theory are cancelled in the cross- 
section of any process if one takes into account the existence in the initial and final states of 
an arbitrary number of the soft quanta as well as the particles having the parallel momenta 
with the same total momentum. The probabilities of these processes integrated over the phase 
space of these additional soft (collinear) quanta in the initial and final states should be added 
to the probability of the initial process. 

As an illustration we consider the model example of the electron-proton (quark) scattering 
and put all the masses equal to zero. We will be interested in the radiative corrections in the 
first order with respect to the strong coupling a^. The corresponding diagrams are shown 
in FigSOl 




a) b) c) d) e) f) 



Figure 40: The process of electron-quark scattering in the first order in a^: ) the Born 
diagram, b)-d) the corrections due to the virtual gluons, e)-f) the corrections due to the real 
gluons 

We have already calculated the matrix elements corresponding to these diagrams, but 
now we proceed in a different way. Since the ultraviolet divergences which appear in the 
diagrams b)-d) are compensated due to the Ward identity in QED {Zi = Z2), all the arising 
divergences are solely infrared and collinear. To extract them we will use the dimensional 
regularization. Then both the divergences are manifested in the form of the poles over e 
and, since we have both of them, there will be poles of the first and the second order. 

We start with the virtual corrections. The diagrams of self-energy c) and d) in the 
massless case are identically zero due to the above-mentioned property of a massless integral 
depending on one argument equal to zero (p^ = on the mass shell). As we explained, here 
one has the cancellation of the UV and the IR divergences. Therefore, all divergences in 
the vertex diagram b) may be considered as infrared. (The UV divergences should cancel 
with the UV ones from the self-energy diagrams and the latter in their turn cancel with the 
IR). The integral for the vertex part is defined by two form factors Fi{q^) and -^2(9^) (19. ip . 
Taking the expression for the vertex function (I3.22p as the starting point, we put m = and 
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go to the mass shell. The result is 
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where instead of the logarithm of the photon mass as the IR regulator we have the pole 
over e. In order to avoid the transcendental numbers, we used the helpful definition of the 
angular measure in the space of 4 — 2s dimensions and multiplied the standard expression by 
r(l — £) / (47r)^. Then the constants like 7^, log{A7r) and ({2) disappear from the intermediate 
expressions. Due to the cancellation of divergences in the final expressions, this redefinition 
does not influence the answer. 

Thus, the cross-section for the diagrams with virtual gluon has the form 
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where the differential cross-section in the Born approximation is given by 
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In the cm. frame s — E^,t — — E^/2(l — cos 6*), it = — E^/2(l + cos 6), where the angle 9 is 
the electron scattering angle. 

Consider now the diagrams with the emission of the real gluons e) and f). Besides the 
squares of each of the diagrams one should also take into account the interference term. The 
calculation in fact repeats that in QED but instead of the photon mass we again use the 
dimensional regularization and do not restrict the integration region over the momentum of 
additional gluon. The calculation is a bit tedious, after contracting all the indices the phase 
integral takes the form 
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It is useful to pass to the spherical coordinates and use the cm. frame. After the integration 
over the phase volume the result can be represented in the form 
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real 
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where the functions /i and /2 in the cm. frame are [x = cos 6) 
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+ ^(1 - x)(3x^ + 15x2 + 77x-31) log(^^) + (l + xf{x^\^x\'i)T{^ 
-12(9a;^ + 2a;+5)Lz2(i^) + ^(l - x)(l + x) (5x^-420; -23) 



(10.9) 



As one can see from the comparison of the cross-sections of the processes with the virtual 
(110.41) and the real gluons (110. 7p . in the sum the second order poles cancel. However, the 
total cancellation of divergences does not happen. The remaining divergences in the form of 
a single pole have a coUinear nature. As was already mentioned, for their cancellation one 
has to define properly the initial states. The point is that the massless quark can emit the 
coUinear gluon which will carry part of the initial momentum and in this case, it is impossible 
to distinguish one particle propagating with the speed of light from the two flying parallel. 

10.2 The quark distributions and the spHtting functions 

To take into account this possibility, let us come back to the scattering process and assume 
that the initial quark has emitted the parallel gluon (see Fig HTj) . The two particles can be 
almost parallel with small relative transverse momentum. The three four-momenta can be 



Figure 41: The diagram corresponding to the splitting of the quark into the quark and the 
gluon 

chosen in the form: 

p=(p;0,0,p), g ^ (zp;p_L,0, zp), /c ^ ((1 - 2;)p; -p_L, 0, (1 - z)p), 

so that all of them obey the condition = = fc^ = with the accuracy up to p^.- 
is helpful, however, to use another method, namely to choose the momenta in such a way 
that they obey the mass shell condition with the accuracy up to p^, but to give up the 
energy conservation in the order of The advantage of this approach consists in the use 
of formulas for the spinors and the polarization vectors on mass shell. Therefore, we choose 
the momenta as follows: 

V = (P;0,0,p), q ^ {zp+ ^;p^,0,zp), k ^ {{l-z)p+ 2(^i^^^p '^ 0, il-z)p). 

The square of the matrix element corresponding to the process of splitting on mass shell in 
this case can be written in the standard form 

|M(g ^ qG)\' = ^^^r(7Wg)Ee*'^e^ (10.10) 

pol 



where the factor 1/2 comes from the averaging over the spin states. Here we must take into 
account the physical polarizations of the gluon only, i.e. 

pol 



which gives 

\M{q^qG)\'=Ag\ 
or, substituting the values of momenta. 
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:io.iii 



|M(g^gG)p = C^^^^^^, z<l. (10.12) 

The obtained expression does not depend on the choice of momenta and has a universal 
character. 

Now one can calculate the cross-section of the process of interest. Graphically, it will 
be the same diagram Fig HQ] ); however, the additional gluon will be referred not to the 
final state but to the initial one. Here we use the standard Feynman rules when the energy 
conservation law is not violated, but the massless particle is slightly off shell. Since in the 
case of interest the quark with momentum q is virtual, it is useful to choose the momenta 
like 

P= (p;0,0,p),g {zp-—^—;p^,0,zp),k^ {{1- z)p +—^—;-p±,0, {l-z)p). 

2[l—z)p 2[l—z)p 

In this case, 

= -T^. (10.13) 

1 — z 

Then the cross-section of the process can be written in the factorized form 

mp) = j^J ^\M,^.on^.n^^)Mpz), (10.14) 

where the factor (^) is due to fact that the cross-section is normalized to the energy of 
initial particles, and we have replaced the quark with the energy p^ by the quark with the 
energy zp^. 

Rewriting the differential d^k in terms of the new variables 

d^k = pdzd'^p± = pdzndp^, 
and substituting the value of the matrix element (110.121) and from (110. 13p . we get 

2 ^ pdzdpl (1 - zf 2p\ 1+^2 

da{p) = Cfz-^ / 7- r 4 jz r- zda{pz) 

167r^ (1 — z)p p\ z[l — z) 1 — z 

as f dzdpll + z^ f,m^\ 
= C,- j ^^daipz). (10.15) 

The integral over the transverse momentum is divergent at zero and this is nothing else but 
the manifestation of the collinear divergence. The upper limit is not of great importance, it 



is restricted by kinematic considerations. We assume that the integration over p\ goes from 
zero to some scale Q^. Later, we will see that one can change this scale analogously to the 
change of the ultraviolet scale /i^. 

To extract the divergence we use the dimensional regularization. Changing the dimension 
of transverse integration from 2 to 2 — 2e one gets 



a 



271 Jo 1 — z Jo p]_ 



e 



At first sight the obtained expression still contains the pole in the integrand as 2; 1. 
However, it only looks like a singularity. It came from the matrix element fll0.12p . which we 
have calculated only for z < 1 and it needs to be redefined for z 1. We will come back to 
this question below and, at first, discuss the interpretation of relation fllO.161) . 

Let us introduce the notion of distribution of the initial quark with respect to the fraction 
of the carried momentum z: q{z). Then the initial distribution corresponds to q{z) = S{l—z), 
and the emission of a gluon leads to the splitting: the quark carries the fraction of momentum 
equal z, while the gluon - {1 — z). The probability of this event is given by the so-called 
splitting functions Pqq{z) and Pqci^ — z). In the lowest order of perturbation theory in as 
the quark and gluon distributions can be written in the form 

g(z,Q^) = + (10-17) 

G{z,Q') = ^\[^^^ P^o{l-zl (10.18) 

where the splitting functions are defined by the corresponding matrix elements one of which 
for Pqq{z) has been calculated in the leading order in a<j earlier (see fll0.12p ). The result has 
the following form: 

P,,{z) = (^^^ + ^5(1 - ^)) , (10.19) 

P,g{z) = + (10.20) 

Note that eq. (110. 19p contains the redefinition of the function Pqq{z) at the point z = 1 
mentioned above, namely the sign " + " should be understood as the following integration 



/o (1 — z)^ Jo (1 — z) 

and the coefficient of the 5-function is defined from the requirement of conservation of the 
number of quarks 

/ q{x,Q^)dz = l I Pqq{z)dz = 0. 

Jo Jo 

Thus, eq. (110. 16p together with the Born diagram can be written as 

da(p) = C dz q(z,Q^) da(pz), (10.21) 
Jo 



where the quark distribution q{z, Q"^) is given by (110.171) . 

It seems strange at first sight that the answer depends on the scale which defines 
the quark distribution. However, it has the physical interpretation. This is the measure of 
collinearity of the emitted gluons that can be distinguished, i.e., it refers to the definition 
of the initial state. In fact, in the massless case one cannot define the initial state that 
contains just the quark, it exists together with the set of collinear gluons. (The same is 
true for the massless electron with collinear photons.) This scale is sometimes called the 
factorization scale, at this scale the scattering cross-section (110. 2ip takes the factorized 
form. The factorization scale can be varied. The dependence of the quark and the gluon 
distributions on the scale is governed by the so-called DGLAP equations well known in QCD. 



10.3 The finite answers 



Thus, besides the two contributions to the cross-section from the virtual and the real gluons 
there is one more contribution related to the splitted initial state (110. 16p . In the lowest order 
of perturbation theory in «<, it can be written as 

where the Born cross-section is given by (110. 5p with the replacement of the initial quark 
momentum p by pz, and the factorization scale is an arbitrary quantity associated with 
the quark distribution function. Note that the scale may depend on z. It is quite natural 
to choose the factorization scale equal to the characteristic scale of the process of interest. 
Thus, in our case this choice corresponds to = —t, where t is the Mandelstam parameter 
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Comparing the obtained expression with (110.40 and (110. 7p we see that the last divergence 
cancels and the final expression for the cross-section of the electron-quark scattering with 
account of possible creation of the gluon in the initial and final states takes the form (x = 



cos( 
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(10.25) 
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This expression is our final answer for tfie cross-section of the physical process of electron- 
quark scattering where the initial and the final state include the soft and collinear gluons. 
It includes also the definition of the initial state and can be recalculated for the alternative 
choice of the factorization scale similar to what happens to the ultraviolet scale which de- 
fines the coupling constant. Thus, we practically deal with the scattering not of individual 
particles but rather with coherent states with a fixed total momentum. Only this process 
has a physical meaning. 

In FigH2| we show the differential cross-section of this process as a function of the 
electron scattering angle: ^^(cos^). We have chosen here the strong coupling as = 0.2, 
and Cf = 4:/3. As one can see, the inclusion of the radiative correction ~ as practically 




Figure 42: The differential cross-section of eg scattering in the Born approximation and with 
allowance for the as correction. On the right plane the same plot is shown in the bigger 
scale 

does not change the result, the difference from the Born approximation is less than a per 
cent, that justifies the use of perturbation theory. 

Let us stress once more that the obtained answer for the cross-section of the observable 
process depends on: a) the ultraviolet subtraction scheme that manifests itself, in particular, 
in the appearance of the ultraviolet scale /x^ (canceled in our case in the lowest order of 
perturbation theory) and b) the definition of the initial coherent state, which manifests itself 
in the appearance of the factorization scale Q'j. The universality in the description of the 
physical processes is based on the fact that choosing the UV and the IR scale one way or 
another and fitting the experimental data of some process, one can then recalculate the 
obtained values of the running coupling and of the quark (lepton) distribution for any other 
choice of the scales. This way the result for the observable quantities does not depend on a 
particular choice of these scales and is universal. 



on 



11 Afterword 



Local quantum field theory, being the mathematical basis of elementary particle physics, is 
the logical continuation of quantum mechanics. It exploits the same basic ideas, but de- 
scribing the system with an infinite number of degrees of freedom permits the creation and 
annihilation of particles in the course of the interaction. The modern formulation is based 
on the interaction representation which assumes the existence of the asymptotic states of 
the free fields. In the S-matric approach we presume that these fields interact in a local 
way in the space-time, and calculating the S-matrix elements one can find the probabilities 
of various processes. The most developed and reliable method of these calculations is the 
perturbation theory in the coupling constant which is similar to the one in quantum mechan- 
ics. However, due to a much more complicated structure of the field theory, the methods 
of perturbation theory encounter problems which have no analogy in quantum mechanics, 
namely the divergence of the appearing integrals for the radiative corrections. We have 
shown in these lectures how one can deal with these divergences which have the ultraviolet 
and the infrared nature and how to get the finite answers for the probabilities of the physi- 
cal processes. We did not aim to prove the main theorems like the Bogoliubov-Parasiuk or 
the Kinoshita-Lee-Nauenberg theorem, but have exemplified how they work. The explicit 
calculations allow one to convince himself in the validity of the final conclusions. 

It should be noted that the formalism of quantum field theory contains the physical prin- 
ciples which we have to follow sometimes not realizing it. Thus, for example, the ultraviolet 
divergences restrict the type of the interaction and, contrary to quantum mechanics, there 
are only a few types of allowed Lagrangians. Not without reason the renormalizability played 
such an important role in the formation of the Standard Model. The other example is the 
notion of the asymptotic states. Even starting with the free fields within the perturbation 
theory, from the requirement of the cancellation of the infrared divergences we come to the 
definition of the physical initial and final states which are essentially the coherent states. 

The very fact that the gravitational interaction does not fit to the general scheme proba- 
bly means that local quantum field theory has a limited applicability and should be replaced 
by a more general construction. It might be nonlocal like in the string theory, or multidi- 
mensional one like in the brane-world theory. However, in any case, in the low energy limit 
one has the local quantum field theory though possibly going beyond the Standard Model 
that we considered here. 
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